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V  ABSTRACT 

^An  Internal  airblast  from  a  missile  explosion  produces  tearing, 
holing,  and  dishing  of  a  ships'  structural  plating. 

In  order  to  analyze  the  residual  strength  of  a  ship  after  an 
Internal  explosion,  a  typical  ship  compartment  was  considered  to  be 
composed  of  flat  square  plates  subjected  to  uniaxial  compression,  and 
the  damage  effects  were  examined  independently.  The  holing  was 
considered  to  be  a  centrally  located  circle,  Plastic  limit  load 
analysis  of  a  flat  plate  indicates  that  ultimate  strength  decreases  in 
proportion  to  increasing  hols  sise.  Elastic-plastic  analysis  indicates 
that  ultimate  strength  is  little  affected  by  hole  size  until  the  holing 
reaches  half  the  plate  width,  and  is  essentially  equal  to  the  ultimate 
load  described  by  the  effective  width  concept.  An  extension  of  the 
concept  of  small  imperfections  to  dished  plates  suggests  that  the 
ultimate  strength  is  reduced  by  roughly  10  percent  from  that  of  a  flat 
plate.  The  extensional  collapse  mode  analysis  describes  folding  and 
strength  beyond  peak  load.  Fully  plastic  crack  propagation  was 
examined,  but  its  limiting  effect  can  not  be  precisely  described. 


A  small  program  of  experimentation  was  developed  to  geln 
qualitative  insights  into  the  interaction  and  cumulative  effect  of  the 
various  damage  modes.  Seven  mild  steel  boxes  were  constructed  of  four 
square  side  plates  with  length- to -thickness  ratios  of  90  and  two 
stiffened  end  plates  to  model  different  damage  effects.  The 
experimental  results  Indicate  that  the  plastic  limit  load  doks^nfiL 
define  the  peak  load  of  slender  plates.  Increased  hole  size  appears't*. 
decrease  ultimate  strength  gradually.  Dishing  on  the  order  of  5  to  15 
plate  thicknesses  reduces  ultimate  strength  by  roughly  lit  £&  percent. 
The  mean  crushing  load  allows  damaged  structural  plating  to  continue  to 
support  loading  on  the  order  of  £  La  IQ  percent  of  the  tensile  yield 
load.  Although  crack  growth  was  not  a  factor  during  testing,  it  could 
be  Important  in  shipboard  structural  plating  which  is  relatively  thinner 
and  harder  than  the  tested  specimens. 
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Most  notation  will  be  specifically  defined  where  it  is  first 
introduced.  The  following  table  will  serve  as  a  reference  guide.  Where 
more  than  one  meaning  has  been  assigned  to  a  symbol,  the  correct  use 
will  be  obvious  from  the  text  or  explicitly  defined. 

A  Cross-sectional  area,  crack  area 

a  Length  of  plate 

b  Width  of  plate 

Effective  width  of  unperforated  plate,  corrected 

effective  width  for  holed  plate 

C  Coefficient  of  deflection  function 

s 

c  Diameter  of  hole 

e  Width  of  tensile  region  in  unperforated  plate 

D  Plate  flexural  rigidity  -  Et*/  [12(1  -  wa>] 
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n  Unit  stress  resultants 
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Perjt  Critical  load  of  unperforated  plate 
P  Compressive  load  removed  by  hole 
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Py  Tensile  yield  load  (squash  load) 

Q,q  Generalized  stress,  strain 

R  Fracture  toughness  (essential  work  of  fracture) 

r  Coordinate,  radius  of  circular  plate  &  tube 
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1 .  INTRODUCTION 


1.1  The  Antlshlp  Missile  Threat 

Modern  naval  warships  oust  be  built  to  survive  and  operate  in  a 
multiple  hostile  threat  environment.  In  addition  to  facing  heavy  seas 
and  storms,  a  warship  may  be  subjected  to  attack  by  a  variety  of  weapons 
such  as  torpedoes,  mines,  gunfire,  bombs,  and  missiles  (Figure  1.1). 
Some  of  the  more  common  and  lethal  weapons  are  antiship  missiles 
equipped  with  semi-armor  piercing  warheads  with  small,  high  explosive 
(HE)  payloads.  Such  antiship  missiles  first  penetrate  a  ship's  outer 
structure  and  then,  by  a  delayed- fuze  mechanism,  detonate  within  the 
confines  of  the  superstructure  or  hull  to  cause  damage  to  the  structure 
and  internal  systems  [1].  The  purpose  of  this  paper  is  to  analyze  the 
structural  ability  of  a  destroyer- sized  warship  to  withstand  the  effect 
of  explosive  loading  associated  with  a  missile  attack. 

The  various  antiship  missiles  produced  by  different  nations  all 
have  similar  components.  A  typical  antlshlp  missile  has  a  propulsion 
soction  (multi-stage  solid  fuel  and/or  a  turbojet  engine  plus  a  booster 
for  surface/submarine  launch),  fuel  tank,  guidance  system,  and  a  warhead 
enclosed  in  an  aerodynamic  cylinder.  Small  antiship  missiles  are  widely 
produced  and  exported  throughout  the  world.  Their  sizes  range  from  the 
now  Infamous  17  foot  French  EXOCET  with  a  200  pound  warhead  down  to  the 
10  foot  Norwegian  PENGUIN  with  a  264  pound  warhead.  Most  antiship 
missiles  attack  their  targets  with  a  sea  skimming  trajectory  which  make 
them  difficult  to  detect  and  track. 

Modern  antiship  missiles'  blast  and  fragmentation  warheads  are 
little  more  sophisticated  than  the  free-fall  bombs  used  against  ships 
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sine*  the  start  of  World  War  Two.  Tha  affact  of  aarlal  bombs  and 
kamikazs  aircraft  attacks  on  destroyers  during  World  War  Two  illustrates 
tha  type  and  extent  of  damage  produced  by  such  weapons.  A  direct  hit  of 
an  aerial  bomb  causes  damage  to  the  above -water  and  underwater  sections 
of  the  hull  (sometimes  with  a  breakoff  of  an  extremity) ,  fragmentation 
damage  to  exposed  stations  and  superstructures,  flooding  of  compartments 
producing  a  change  in  draft,  and  tha  outbreak  of  fires.  The  main 
propulsion  spaces,  electrical  equipment,  and  the  ship  guns  could  be  put 
out  of  action  as  a  result  of  a  high  explosive's  blast  and  fragmentation 
damage.  Aerial  bomb  explosions  put  destroyers  out  of  action  basically 
because  of  damage  to  guns  and  technical  equipment  (about  50  percent  of 
the  cases,  of  which  almost  25  percent  were  due  to  shock  vibrations)  and 
because  of  loss  of  watertight  Integrity  (35  percent) .  Fires  and 
explosions  and  destruction  of  the  overall  strength  of  the  hull  of  the 
ship  had  less  significance.  The  explosive  action  of  kamikaze  aircraft 
most  often  caused  destroyers  to  lose  mission  capability  rather  than  to 
sink  [2].  Of  30  U.S.  destroyers  sunk  during  World  War  Two  due  to 
gunfire,  bomb  or  kamikaze  attack,  only  six  resulted  from  loss  of 
structural  integrity  of  the  hull  (Figures  1.2,  1.3)  [3]. 

Significant  ship  structural  damage  due  to  weapons  effects  has  been 
scant  since  World  War  Two  and  has  generally  not  been  as  well  documented 
[3].  During  the  recent  Falklands  Conflict,  however,  the  British  Royal 
Navy  faced  a  multiple  hostile  threat  environment.  HMS  Ardent .  a  3250 
ton  frigate  suffered  four  aerial  bomb  hits  before  it  was  "mission 
killed"  and  required  three  more  hits  before  it  sank.  HMS  Antelope, 
sister  ship  of  HMS  Ardent .  sank  after  one  500  lb.  bomb  exploded  and 
started  an  uncontrollable  fire  that  caused  several  secondary  explosions. 


13 


HMS  Sheffield  vn  hit  by  an  EXOCET  antiship  missile  which  failed  to 
explode,  but  which  penetrated  and  ignited  the  ship's  main  fuel  supply 
tank  with  its  still  burning  rocket  fuel.  The  ensuing  fire  and  smoke  was 
uncontrollable.  The  15,000  ton  roll-on/roll-off  container  ship  Atlantic 
Conveyor  was  also  lost  due  to  a  single  EXOCET  missile  hit  which  produced 
uncontrollable  fires  [4] . 

1.2  Modem  U.S.  New  Warship  Design 

United  States  warship  structures  have  always  been  implicitly 
designed  with  survivability  in  mind.  Incorporating  lessons  learned  from 
World  War  Two  into  such  design  practice  as  crack  arrestors,  requiring 
gradual  transitions,  using  continuous  welds,  providing  hull  girders  with 
a  1.0  ton  per  square  inch  reserve  primary  strength  and  limiting  design 
primary  stress  throughout  the  structure,  has  produced  more  survlvable 
warship  designs.  David  W.  Taylor  Research  and  Developmenc  Center  has 
been  developing  material  fracture  control  technologv.  Strict 
fabrication  requirements  minimize  the  deleterious  effects  of  poor 
structural  detailing,  failed  welds,  ctress  concentrations,  etc.  in  high 
Intensity  dynamic  pressure  environments  and  also  contribute  to  the 
survivability  of  the  warship  [3] .  Current  techniques  for  predicting 
weapon  effects  on  structures  are,  however,  limited.  To  simplify  the 
determination  of  residual  strength  after  an  attack,  the  Naval  Sea 
Systems  Command  considers  all  damaged  structure  to  be  "ineffective"  [3]. 

Small,  unarmorod,  warships  such  as  destroyers  are  likely  to  suffer 
local  structural  failure  from  the  initial  explosion  and  be  faced  with 
spreading  fire  causing  further  damage  and  secondary  explosions.  The 
large  amount  of  flammable  material  on  modern  warships  make  them  more 
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susceptible  to  firs  damage  than  World  War  Two  warships.  Ona  of  the 
lessons  learned  from  the  Falklands  Conflict,  and  further  emphasized  in 
the  Persian  Gulf,  was  that  aluminum  structures  failed  because  of  the 
intense  heat  produced  by  fires.  Although  steel  structures  may  be 
weakened  or  distorted  by  shipboard  fires  on  the  order  of  1800°F  to 
2000°F,  aluminum  structures  will  rapidly  sag  and  melt  at  such 
temperatures  [5].  Since  World  War  Two,  small  warships  were  built  with 
aluminum  superstructures  to  reduce  topside  weight  and  associated  costs. 
The  U.  S.  Navy's  new  DDG-S1  Arlelgh  Burke  Class  destroyer  incorporates 
improved  design  practices,  all  steel  construction,  and  damage  control 
features  such  as  shock  mounted  and  blest  hardened  structures  and 
equipment  along  with  redundant  firefighting  capabilities.  Such  warships 
can  undoubtably  survive  the  secondary  effects  of  missile  damage  better 
than  the  the  destroyers  of  World  War  Two  or  the  British  warships  in  the 
Falklands  Conflict. 

Contemporary  tactical  doctrine  calls  for  coordinated  saturation 
attacks  by  antiship  missiles  to  overcome  a  warship's  active  defenses  and 
ensure  several  missile  hits  14].  Even  with  severe  battle  damage,  e 
hardened  structure  provides  a  warship  with  a  "keep  afloat"  capability. 
To  Improve  this  passive  defense  capability  against  the  growing  antiship 
missile  threat  it  is  important  to  understand  how  a  warship's  structural 
plating  is  affected  by  the  damage  produced  by  an  internal  explosion. 
The  capacity  of  a  warship's  structure  to  absorb  and  contain  missile  hits 
Is  vital  to  its  ability  to  survive  and  operate  in  a  hostile  environment. 

1.3  Local  Pamagc  Effect!  And  M&fihad  o£  Analyli 

An  internal  airblast  from  a  missile  explosion  imparts  reverberating 


high- intensity  pressure  pulsaa  against  ship  structural  plating.  The  HE 
detonation  produces  gasses  with  a  temperature  on  the  order  of  5500°F  and 
pressure  reaching  750,000  psi.  The  total  impulse  is  a  function  of 
missile  momentum  and  explosive  energy.  After  passage  and  decay  of  the 
initial  shock  wave,  confinement  of  the  explosion  gases  leads  to  a 
pressure  loading  with  a  typical  duration  on  the  order  of  several 
seconds  [6]. 

Structural  damage  is  caused  by  a  combination  of  (1)  holing  of 
Internal  structural  plating,  (2)  blast  overpressure  within  confined 
spaces,  (3)  fragments  from  warhead  break-up,  (4)  fragmentation/spalling 
of  interior  bulkheads  and  decks,  and  (5)  fire  [1].  These  modes  of 
failure  interact  and  may  cause  structural  collapse  of  a  local  or  global 
character.  The  effect  of  fragmentation  on  structural  integrity  is 
generally  insignificant  and  is  not  evaluated.  The  subsequent  effect  of 
heat  and  fires  produced  by  missile  damage  on  structural  integrity  has 
been  analyzed  by  Hlddeman  [7].  The  structural  effects  of  supporting 
stiffeners  (girders,  beams,  etc.)  have  been  examined  by  Faulkner  [8]  and 
others.  This  paper  is  specifically  concerned  with  the  residual  strength 
in  a  ship's  structural  plating  which  has  suffered  an  internal  airblast. 

The  dynamic  interaction  of  damage  effects  within  a  complex  ship 
structure  is  not  well  understood.  In  order  to  analyze  this  situation,  a 
typical  warship  compartment  structure  is  considered  to  be  composed  of 
flat  square  plates  subjected  to  quasi-static  uniaxial  loading,  and  the 
damage  effects  are  examined  independently.  Structural  damage  consists 
of  some  combination  of  tearing,  holing,  and  deflection. 

The  extent  of  holing  caused  by  the  missile  explosion  is  based  upon 
U.S.  Navy  experiences  [1],  Holing  is  caused  by  the  high 
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velocity  gradients  across  the  plate  surface  (due  to  localized  pressure 
loading)  [6].  Roberts  and  Azlzlan  [9)  have  shown  analytically  and 
experimentally  that  hole  shape  has  little  effect  on  structural  strength. 
As  a  reasonable  simplification,  therefore,  the  holing  is  considered  to 
be  a  centrally  located  circle. 

As  an  initial  approximation  of  residual  strength  after  damage,  a 
plastic  limit  load  analysis  and  elastic-plastic  buckling  analysis  are 
developed  for  a  flat  plate  with  a  circular  hole.  The  dishing  of  a  flat 
plata  under  impulsive  loading  is  examined  to  determine  the  relationship 
between  lateral  deflection  and  blast  overpressure.  The  residual 
strength  of  a  dished  plate  is  examined  by  an  extension  of  the  concept  of 
initial  imperfections  and  by  the  recently  developed  extenslonal  collapse 
mode  analysis.  The  limiting  effect  of  mixed  mode  crack  propagation  upon 
residual  strength  is  then  examined.  This  progressive  analytic  approach 
is  presented  in  Figure  1.4. 

A  modest  program  of  experimentation  is  conducted  to  determine  the 
relationship  between  deflection,  holing,  crack  propagation  and  rasidual 
strength.  An  analysis  is  then  made  of  the  residual  strength  in  a  modern 
U.S.  Navy  destroyer  (DDG-51  Arleigh  Burke  Class)  hit  by  one  antiship 
missile  amidships. 
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Figure  1.1: 


Categorization  of  naval  antiahip  weapons  [3] 
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Figure  1.2:  Diagram  of  battle  damage  to  USS  Brownson  (DD-518)  which 
sank  due  to  loss  of  structural  integrity  [3] 
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Figure  1.3:  Summery  of  bottle  damage  to  several  World  War  Two  U.  S. 
destroyers  which  sank  due  to  loss  of  structural  integrity  [3] 


Figure  1.4:  Analysis  of  battle  damage,  (a)  Real  damage  (b)  Uniaxial 
compressive  loading  of  a  flat  plate  with  a  central  circular  hole  (c) 
Uniaxial  compressive  loading  of  a  dished  plate  (d)  Uniaxial  compressive 
loading  of  a  dished  plate  with  central  circular  hole  and  initial  cracks 
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2.  PLASTIC  LIMIT  LOAD  OF  A  FLAT  PLATE 

2.1  Introduction 

Shipboard  structural  members  such  as  bulkheads  and  decks  can  be 
modeled  as  thin  plane  sheets  of  homogeneous  material.  Such  a  structural 
member  la  often  acted  upon  by  forces  In  Its  plane  which  are  much  larger 
than  the  normal  and  shearing  forces  acting  on  Its  plane.  If  the  sheet 
lies  In  the  x-y  plane,  the  non-vanishing  stress  and  strain  components 
are  a  ,  a  ,  a  and  <  ,  <  ,  <  respectively.  The  loads  which  causa 

xx  yy  xy  xx  yy  xy  y 

plastic  collapse  (yield  failure)  In  a  structure  can  be  analytically 
described  through  the  use  of  lower  and  upper  bound  theorems. 

Structural  damage  consists  of  some  combination  of  tearing,  holing, 
and  deflection.  The  compliance  of  a  ship  generally  imposes  a  uniform 
displacement  (which  some  authors  describe  as  constant  strain)  boundary 
condition  on  the  loaded  edges.  A  rough  upper  bound  limit  load  analysis 
of  the  effects  of  fully  plastic  fracture  in  a  dished  plate  with  a 
central  circular  hole  is  presented  in  Chapter  6.  As  an  initial  analysis 
of  the  effects  of  damage,  however,  a  plastic  limit  load  analysis  of  a 
flat  plate  with  a  central  circular  hole  subjected  to  uniaxial 
compressive  loading  is  presented  here. 

2.2  Physical  Basis  q£  Plasticity 

Plastic  deformation  occurs  when  a  structure  experiences  loads  which 
produce  stresses  above  the  yield  strength.  The  stress -strain  curve  for 
mild  steel,  shown  in  Figure  2.1,  illustrates  actual  versus  perfectly 
plastic  deformation.  Microscopically,  plastic  deformation  consists 
primarily  of  shear  relative  to  crystallographic  slip-planes,  with  no 
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volume  change.  Three  consequence*  of  the  micromechanics  carry  through 
to  macromechanics:  (1)  incompressibility,  (2)  current  stress  is  related 
to  the  increment  of  strain,  not  to  the  total  strain  as  is  the  case  with 
elasticity,  and  (3)  the  associated  flow  rule.  A  continuum  theory  of 
plasticity  can  only  be  applied  to  relatively  large  regions  as  compared 
with  the  mean  free  path  of  dislocations  and  the  grain  size  (0.001  to 
0.1mm),  or  perhaps  to  coherent  slip  bands. 

There  is  a  relation  between  the  yield  locus,  which  determines  the 
state  of  stress  at  which  plastic  flow  can  occur,  and  the  equations 
relating  stress  and  strain  increment  for  the  subsequent  plastic  flow. 
The  associated  flow  law  establishes  a  one-to-one  correspondence  between 
the  state  of  stress  and  strain  rate  if  the  yield  locus  is  smooth  with  a 
finite  curvature  at  each  point.  It  has  been  shown  to  hold  not  only  at 
the  microscopio  level  but  for  forces  and  momenta  on  structures  as  a 
whole.  Once  the  yield  strength  has  bean  reached,  the  current  state  of 
stress  determines  the  ratios  between  the  strain  Increments  but  does  not 
define  the  strain.  A  hole  in  a  tensile  structure  under  plastic  loading 
will,  for  non-hardening  plasticity,  localize  the  strain  in  a  narrow  band 
and  that  band  will  extend  all  the  way  to  the  nearest  boundary  of  the 
part  [10], 

2.3  Limit  Laid  Analyili 

As  a  structure  experiences  general  plastic  deformation,  its 
load-deflection  curve  becomes  nearly  horizontal.  The  limit  load  Is  the 
horizontal  limit  approached  by  the  load- deformation  curve.  Although  it 
can  be  very  difficult  to  analytically  determine  the  Unit  load  of  a 
structure,  it  is  comparatively  simple  to  determine  upper  and  lower 
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bound*  to  th*  limit  load. 


According  to  tha  lover  bound  the  or  an,  any  load  which  ia  in 
equilibrium  with  a  stres*  field  *uch  that  the  yield  condition  1*  not 
violated  and  the  equilibrium  equation*  are  satiafied  is  a  lover  bound  on 
the  collapse  load.  The  equations  of  equilibrium  for  plane  stress  are: 

6a  da  do  da 

- ”  +  _i!  _  0  — =1  +  - 2i-0  .  (2.1) 

fix  fly  fix  fly 

There  are  several  way*  to  express  the  yield  condition.  Hodge  [11]  used 
the  Tresca  yield  criterion  which  can  be  expressed  as: 
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and  Gordon  and  McCrum  [12}  performed  a  limit  load  analysis  using  th* 
Mises  yield  criterion  which  can  be  expressed  as: 
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where  cry  is  th*  material  yield  strength.  In  constructing  lower  bounds, 
it  is  frequently  useful  to  consider  discontinuous  stress  fields  across  a 
plan*  with  normal  n  and  in-plan*  coordinates  s,t.  Then  th*  components 
acting  across  th*  surface  a  ,a  ,a  must  be  continuous  for 
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equilibrium,  while  th*  ln-plan*  components  a  ,a  ,a  need  not  be. 

ti  t  IB  B  t» 

An  upper  bound  on  th*  limit  load  can  be  found  by  considering  any 
incompressible  velocity  field  which  satisfies  all  velocity  boundary 
conditions.  It  is  frequently  useful  to  consider  discontinuous  velocity 
fields  when  evaluating  upper  bounds,  which  is  perminslbl*  as  long  as  th* 


work  associatad  with  such  •  discontinuity  la  accounted  for.  The  yield 
locua  is  a  contour  of  constant  value  of  the  flow  function  which 
characterizes  the  yield  condition.  Analysis  of  the  yield  hexagon 
reveals  that  the  internal  energy  is  the  product  of  the  yield  stress  by 
the  absolutely  largest  of  the  three  principal  strain  rates .  Any  load 
which  produces  an  external  energy  rate  in  excess  of  this  value  is  an 
upper  bound  on  the  Unit  load.  For  the  case  cf  uniform  traction  (which 
some  authors  describe  as  constant  stress),  the  upper  bound  to  the  limit 
load  F+  can  be  expressed  as: 


p+  _  6U* 
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where  A  is  the  gross  cross-sectional  area,  SWP  la  the  increment  of 
plastic  work,  and  Sup  ia  the  increment  of  deflection  along  the  line  of 
loading.  For  the  case  of  uniform  displacement,  the  upper  bound  to  the 
limit  load  can  be  simplified  ae: 
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2.4  Piastre  Limit  Load  of  g  Holed  Plate 

The  limit  load  of  a  flat  square  plate  with  a  central  circular  hole, 
subjected  to  uniaxial  compressive  loading,  can  be  determined  within 
reasonable  upper  and  lower  bounds.  Hodge  [11]  analyzed  the  tensile 
limit  load  of  a  flat  square  plate  with  a  central  circular  hole,  simply 
supported  edges,  and  uniform  traction  boundary  conditions,  Figures 
2.2a,b,c  demonstrate  that  similar  but  orthogonal  compressive  or  tensile 
loads  will  produce  the  same  rigid-plastic  failure.  This  can  also  be 
demonstrated  mathematically  by  changing  the  signs  of  loading  and 
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Associated  deflection  in  Eqs.  2.2  end  2.3.  The  rets  of  external  work 
required  to  produce  similar  deflections  is  the  same  for  both  compressive 
and  tensile  loads. 

Hodge  constructed  a  kinematically  admissible  velocity  field  by 
considering  a  sliding  of  the  slab  out  of  its  plane  along  a  plane  through 
the  axis  of  symmetry.  From  the  Tresca  yield  criterion,  an  upper  bound 
P*on  the  limit  load  of  a  flat  plate  of  width  b  and  with  a  central 
circular  hole  of  diameter  c  is: 

P+  -  pY  (1  *  |  )  (2.4) 

where  is  the  tensile  yield  load  defined  as; 

Py  “  *y  A  .  (2.3) 

For  reasonably  large  values  of  ^  ,  however,  a  better  upper  bound  value 
must  be  based  upon  the  bending  model,  utilising  the  oonoept  of  a  plastic 
hinge,  If  £  is  almost  1,  then  the  holed  plate  could  easily  collapse  by 
means  of  yield  hinges  at  its  four  narrowest  points.  The  yield-hinge 
concept  can  be  used  as  a  guide  to  obtain  a  reasonable  velocity  field. 
Tho  upper  bound  limit  load  for  ^  ^  <  1  thus  becomes: 

»*-/»•*!♦  ’'I’1  •  I  •  «•«> 

A  lower  bound  for  this  problem  was  obtained  by  considering  the 
discontinuous  stress  field  shown  in  Figure  2.3.  The  lower  bound  limit 
load  for  0  <  |  <  0.443  ii: 

P- .  b  <2.7) 

C 

where  (  is  a  distance  parameter  chosen  to  reflect  a  discontinuous  stress 
field  for  maximuum  F*.  For  larger  size  holes,  with  0.443  <  g  <  1  ,  the 
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lower  bound  limit  load  it: 

p-.bii_^  .  (2.., 

Tha  upper  and  lower  bounds  are  combined  in  Figure  2.4  to  describe 
the  plate  failure  loading  for  different  size  holes  with  uniform 
traction.  Tha  upper  and  lower  bounds  derived  by  Gaydon  and  McCrum  [12] 
using  the  Mises  criterion,  presented  in  Figure  2.5,  are  very  similar  to 
those  derived  by  Hodge.  Gaydon  [13]  also  performed  a  biaxial  limit  load 
analysis  of  this  case  with  roughly  similar  results. 

For  uniform  displacement  of  the  loaded  edges,  which  more  closely 
approximates  shipboard  structural  boundary  conditions,  limit  load 
analysis  with  Tresca  criterion  turns  out  to  provide  an  exact  solution. 
The  upper  bound  displacement  field  is  described  by  two  45  degree 
seotions  across  the  mid-plane,  as  shown  in  Figure  2.6.  The  lower  bound 
is  described  by  the  constant  stress  regions  outside  the  central  hole,  as 
also  shown  in  Figure  2.6.  The  lower  and  upper  bounds  coincide  to  define 
the  limit  load: 

P  -  P,  [l  -  f)  .  (2.9) 

This  simple  linear  relationship,  which  corresponds  to  the  upper  bound 
for  small  holes,  is  depicted  on  Figure  2.4  along  with  the  uniform 
traction  upper  and  lower  bounds.  Figure  2.4  suggests  that  shipboard 
compliance,  which  restricts  non-uniform  displacement  of  structural 
plating,  improves  the  strength  characteristics  of  the  plating. 

A  separate  method  of  determining  collapse  load  consists  of  simply 
evaluating  the  stress  experienced  by  the  net  plate  structure  at  the 
mid-plane  and  comparing  it  to  the  yield  strength  for  the  material. 
Stress  concentration  effects,  in  general,  apply  only  to  elastic  behavior 
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and  ara  not  consldarad  In  thia  approach.  In  tha  case  of  a  flat  plate 
with  a  circular  hole,  the  stress  experienced  in  the  remaining  structure 
is  increased  and  the  uniaxial  load  required  for  yield  is  correspondingly 
decreased  as: 

f  -  P,  [l  -  |)  ■ 

This  simple  "engineering"  approximation  is  identical  to  the  uniform 
displacement  limit  load. 

The  limit  load  theorems  are  based  on  current  geometry.  Structural 
deformations  may  significantly  Influence  the  boundary  conditions  and 
equations  of  equilibrium,  producing  discrepancies  between  the  actual 
limit  load  and  that  based  on  the  original  shape.  Ductile  structural 
materials  may  not  behave  in  a  perfectly  plastic  manner.  Thin  plates,  in 
particular,  can  experience  considerable  deformation  in  the  elastic  range 
which  distorts  tha  structure  from  the  assumed  plane  shape  prior  to 
plastic  collapse  [14].  Despite  these  limitations,  if  elastic-plastic 
buckling  has  not  occurred,  plastic  limit  load  analysis  provides  a 
reasonable  first  look  at  tha  effects  of  holing. 
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Figure  2.1:  Greph  of  typical  etreis  versus  strain  for  mild  steal  [11] 
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Figure  2.2:  Square  plate  with  a  circular  hole,  (a)  Tensile  load  (b) 
Compressive  load  (c)  Uniform  traction  hinge  type  mode  of  collapse  [11] 
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Figure  2.4:  Graph  of  bounda  for  a  square  plate  with  a  central  circular 
hole  using  Trssca  criterion  (11} 
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Figure  2.5:  Graph  of  uniform  traction  bounds  for  a  square  plate  with  a 
central  circular  hole  using  Misss  criterion  [12] 


Figure  2.6:  Stress  (a)  and  displacement  (b)  fields  for  a  square  plate 
with  a  central  circular  hole  (uniform  displacement) 


3.  ELASTIC - FLASTI C  BUCKLING  AND  POST-BUCKLING  ANALYSIS 


3.1  intrgdufiiign 

Elastic-plastic  buckling  is  a  typical  failure  pattern  of  thin  flat 
plates  under  in-plane  compressive  loading.  The  ultimate  strength  of 
structural  members  defines  their  ability  to  continue  providing  support 
after  suffering  damage,  and  is  reached  after  compressive  loading  has 
caused  stresses  above  the  yield  strength.  The  effective  width  concept 
for  unperforated  flat  plates  describes  the  limiting  effect  of 
slenderness  upon  ultimate  strength.  The  buckling  of  flat  plates 
weakened  by  circular  holes  has  been  fairly  extensively  studied 
theoretically  and  experimentally.  Kurshln  and  Matveev  [15]  made  a 
comparative  analysis  of  five  theoretical  and  three  experimental  studies 
of  the  stability  of  a  square  plate  with  a  central  circular  hole.  The 
best  agreement  between  theoretical  analysis  and  experimental  results  was 
achieved  by  Ritchie  and  Rhodes  [16],  Ritchie  and  Rhodes  examined  square 
and  rectangular  plates,  with  simply  supported  unloaded  edges,  under 
uniform  displacement  compressive  loading  for  both  buckling  and 
post-buckling  conditions.  The  effective  width  concept  discussed  below 
supports  their  conclusion  that  post-buckling  strength  is  not 
significantly  affected  by  small  or  medium  sized  holes. 

3.2  EffBCtlVB  Width  Concept 

For  a  plate  subjected  to  uniaxial  compressive  loading,  in  general, 
only  the  parts  of  the  plate  close  to  the  unloaded  edges  experience 
plastic  yield.  The  stress  distribution  shown  in  Figure  3.1  is  due  to 
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the  "bowing"  effect  in  the  central  son*  of  the  plete.  In  1932,  von 
KArmAn  postulated  that  the  maximum  load  a  plate  with  simply  supported 
edges  can  sustain  is  reachad  when  the  edge  stress  equals  the  yield 
strength  [17].  The  effective  width  of  the  plate  b#,  also  shown  in 
Figure  3.1,  can  be  described  as: 
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where  Ap  is  the  plate  reference  slenderness: 
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and  where  h  is  the  plate  thickness,  and  K  is  the  plate  buckling 
coefficient: 


K  - 


(3.3) 


where  a  and  b  are  the  plate  dimensions,  and  m  is  a  wavelength  parameter 
chosen  to  yield  the  smallest  eignenvalue.  A  modern  warship's  strucural 
plating,  disregarding  the  effect  of  stiffeners,  has  a  typical 
length- to -thickness  ratio  (b/h)  on  the  order  of  200  to  300.  Vlth  a 
reference  slenderness  ranging  from  4  to  8,  such  plating  is  relatively 
very  thin.  By  combining  Eqs.  3.1  and  3.2,  the  effective  width  can  be 
described  as: 


(3.4) 


With  simplified  dimensions  of  8  feet  by  8  feet,  and  a  thickness  of  half 
an  inch,  warship  structural  plating  has  an  effective  width  of  roughly 
1.6  feet.  For  plates  with  high  slenderness  values,  the  effective  width 
is  small  and  damage  in  the  central  zone  has  little  impact  on  residual 
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strength. 


3.3  Buckling  and  Post-Buckling  Analysis 

It  is  possibls  to  analyze  ths  initial  instability  of  a  flat  plats 
with  a  central  circular  hole,  according  to  Ritchie  and  Rhodes  [16],  as 
simply  that  of  a  flat  plate.  The  strain  energy  U  stored  In  a  plate  can 
be  written  as  the  sun  of  the  bending  and  membrane  strain  energies: 


U  -  11  +  U  (3.5) 

b  ■ 

where  U  is  the  bending  strain  energy  and  U  is  the  membrane  strain 
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energy.  These  strain  energies  can  be  expressed  in  terms  of  the 
out-of-plane  deflection  w  and  the  membrane  stress  resultants  and 

N  as: 
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where  D  is  the  bending  rigidity  of  the  plate,  w  is  out-of -plane 

displacement,  and  v  is  Poisson's  ratio.  Ritchie  and  Rhodes  obtained 

the  stress  resultants  N  ,  N  and  N  in  each  element  with  respect  to  a 

xx  yy  xy 

unit  applied  displacement  of  the  loaded  edges.  The  magnitude  of  all 
stress  resultants  were  then  related  to  the  loaded- edge  displacement  by  a 
magnification  factor,  and  the  strains  related  to  the  stresses  by  Hooke's 
law.  The  deflections  wars  represented  by  a  single  sine  wave  in  the 
loaded  direction,  which  is  a  reasonable  approximation  for  plates  with 
small  holes  whose  behavior  is  similar  to  that  of  unperforated  plates, 
but  may  be  less  accurate  for  plates  with  large  holes. 
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Ritchie  and  Rhodes  rewrote  ths  strain  energy  in  matrix  form  in 
order  to  solve  for  the  buckling  modes  and  loads.  The  non* trivial 
solutions  of  the  simultaneous  equations  minimising  the  strain  energy 
were  the  eigenvalues  of  the  edge  displacement  coefficient  4  end  the 
corresponding  eigenvectors.  The  solutions  also  satisfy  the  criterion 
for  neutral  equilibrium.  Thus  the  smallest  eigenvalue  and  its 
corresponding  eigenvector  give  the  buckling  load  or  displacement  and  the 
buckling  mode.  Ritchie  and  Rhodes  used  a  finite  element  analysis,  with 
a  mesh  of  97  triangular  constant-strain  elements  over  one  quarter  of  the 
plate,  to  obtain  the  stress  distribution  for  both  uniform  edge 
displacement  and  uniform  edge  stress  cases.  The  uniform  displacement 
and  uniform  traction  loading  curves  thus  obtained  are  presented  in 
Figure  3.2. 

The  post-buckling  behavior  was  examined  by  a  oomblned  semi-energy 
and  finite -element  method.  Zt  was  assumed  that  post-buckling 
out-of-plane  deflections  remain  of  the  same  form,  but  with  varying 
deflection  magnitude,  as  those  obtained  from  the  buckling  solution. 
This  assumption  was  considered  reasonable  for  plates  without  holes  if 
the  applied  load  was  less  than  about  twloe  the  buckling  load.  The 
presence  of  local  buckles,  however,  causes  changes  in  the  stress  system 
throughout  the  plate.  These  changes  can  be  evaluated  from  the  solution 
of  von  Kdrmdn’s  compatability  equation,  which  relates  ln-plane  and 
out-of -plane  deformations  in  a  plate  by  means  of  a  stress  function  F,  as 
follows : 
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A  second  system  of  stresses  was  darivad  by  solving  von  Kdrmdn's 
compatability  aquation,  and  applying  tha  boundary  conditions  at  tha 
plata  supportad  adgas.  A  third  strass  systam  was  obtalnad,  using  flnlta 
elements,  by  hypothetically  applying  a  loading  system  around  the  hole 
edge  to  counteract  tha  stresses  produced  by  tha  second  systam.  Tha 
complete  membrane  strass  systam  in  tha  plata  was  obtalnad  from  tha 
superposition  of  these  three  systems: 


+  c,  C  +  ) 


yy 


"  ^  n  +  Ca  f  n  ,  +  n  ) 

t  »  yyl  yy3  J 


yyi 


n  -  ^  n  .  +  Ca  f  n  ,+n  .  1 
xy  xyl  *  w  *yS  xy*  * 


(3.8a) 

(3.8b) 

(3.8c) 


lyi  xyi 

where  tha  subscripts  i,  a,  and  a  refer  to  the  respeotlve  stress  system. 

The  principle  of  minimum  potential  energy  was  applied  by  Ritchie 
and  Rhodes  to  their  post-buckling  analysis  in  order  to  determine  the 
buckled  form.  A  change  in  natural  wavelength  of  tha  rectangular  plate 
buckled  form,  from  two  to  three  half ‘wavelengths,  occurs  if  the  hole 
diameter  is  greater  than  0.35  times  the  plata  width.  Brush  and  Almroth 
[18]  developed  a  plot  of  buckling  coefficient  K  versus  aspect  ratio 
(a/b) ,  as  shown  in  Figure  3.3,  which  presents  eigenvalue  curves  for 
critical  values  of  axial  stress  for  simply  supported  unperforated  plates 
subjected  to  ln*plane  compressive  loading.  A  simply  supported  square 
plate  would  be  expected  to  have  a  buckled  form  of  one  half -wavelength. 


3.4  UltllMt*  Strength  A  Holed  Plate 

A  structure  can  continue  to  support  loading  after  buckling.  As 
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with  von  KArain' a  analysis,  tha  failura  of  a  plats  undar  a  compressive 
sdgs  load  is  assumed  to  occur  whsn  tha  maximum  strass  in  tha  loadad 
direction  raachss  tha  yield  stragth.  For  a  simpla  analysis,  Ritchie  and 
Rhodes  defined  the  ultimate  load  as: 


P  ,  -  *w  b  h 

u  It  Y  • 


(3.9) 


where  is  the  effective  width  corresponding  to  tha  yield  stress  oy, 
which  is  similar  to  Eq.  (2.5).  The  effect  of  a  hole  in  a  buckled  plate 
on  the  residual  strength  can  be  approximated  by  redistributing  tha 
compressive  stresses  released  by  tha  introduction  of  tha  hole  over  tha 
remaining  region.  This  produces  an  increase  in  maximum  strass  for  a 
given  load,  and  a  corresponding  reduction  in  affective  width.  The 
post-buckling  stress  distribution  can  be  represented  by  a  simple  sine 
wave  as  shown  in  Figure  3.4  (which  is  similar  to  that  shown  in  Figure 
3.2).  The  increase  in  maximum  stress  due  to  the  hole  is  calculated  by 


redistributing  the  compressive  stresses  "removed"  by  the  hole: 
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where  c  is  tha  region  carrying  tensile  stresses  in  the  unperforated 

plate,  c  is  the  diameter  of  the  hole,  and  P  is  the  compressive  load 

n 

removed  by  the  hole.  The  change  in  stress  is: 

P. 
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For  a  given  external  load  P: 
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where  (a  +  L  a  )  is  the  maximum  stress  after  the  introduction  of  the 
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hole  and  b  is  the  corrected  effective  width.  The  corrected  effective 


width  can  than  ba  written  at: 


b 


Figure  3.5  was  produced  using  the  relationship  between  applied  load  and 
affective  width  separately  derived  by  Rhodes  and  Harvey  [19].  The 
ultimate  load  for  plates  with  various  sire  holes  was  then  calculated 
and  plotted  in  Figure  3.6,  which  indicates  that  ultimate  strength  is 
essentially  unaffected  by  the  presence  of  small  holes  (diameter  less 
than  half  the  width)  and  is  roughly  equal  to  three  times  the  critical 
buckling  load  P  . 

No  attempt  was  made  to  redistribute  tension  stresses  released  by 
the  introduction  of  the  hole  into  the  analysis.  This  is  a  conservative 
approximation  because  the  effect  of  redistributing  high  tensile  stresses 
could  indicate  that  the  effective  width  of  a  plate  with  holing  is 
greater  than  that  of  an  unperforated  plate,  While  it  la  possible  that 
this  could  Indeed  be  the  case  in  practice,  Ritchie  and  Rhodes  suggest 
that  it  would  be  dangerous  to  advocate  such  a  hypothesis  for  use  in 
design  without  experimental  proof. 

Roberts  and  As is lan  [9]  have  examined  this  question  further  and 
state  that  for  uniaxial  and  biaxial  compression  the  buckling  load  may 
increase  with  increasing  size  of  the  central  circular  hole  due  to  the 
redistribution  of  membrane  stresses  towards  the  edges  of  the  plate.  The 
Increase  in  buckling  load  is  dependent  upon  the  boundary  conditions  and 
appears  significant  only  for  plates  with  clamped  edges.  Roberts  and 
Azizian  conclude  that  "even  though  the  buckling  load  may  Increase  with 
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increasing  size  of  the  hola,  It  la  not  expected  that  thara  will  ba  a 
corresponding  lncraaaa  in  tha  ultimata  or  oollapaa  load." 


3.5  Comparison  With  Experimental  Raaulta 

Tha  experimental  invaatigation  of  Ritchia  and  Rhodaa  [16]  consisted 
of  a  sariaa  of  taata  on  aquara  and  raetangular  hot-dip  galvanlzad  mild 
atael  platas  to  collact  information  on  buckling  loads,  out-of-plana 
daf lections,  oollapsa  loads,  and  attain  distributions  along  tha  axaa  of 
the  plates.  Plata  dimensions  ware  nominally  254mm  wide  by  254mm  or 
508mm  long  (square  and  rectangular  respectively)  and  1.57mm  thick.  Tha 
unloaded  vertical  edges  of  the  plate  were  supported  in  tha  frame  by 
knife  edges  and  loaded  edgaa  were  supported  in  rollers  resting  on  needle 
bearings  contained  in  the  blocks ,  through  which  the  load  was 
transmitted.  Tha  specimens  ware  thus  simply  supported  on  all  edges  and 
subjected  to  uniform  displacement  loading.  The  buokllng  loads  were 
estimated  by  obtaining  the  intersection  of  the  prc-  and  post-buckling 
slopes  of  the  load  membrane  strain  plots.  All  plate  specimens  were 
loaded  under  uniform  displacement  conditions. 

Ritchie  and  Rhodes  normalized  their  analytic  solutions  and 
experimental  results  with  respect  to  the  critical  buckling  load.  From 
Timoshenko  and  Gera  [20],  the  critical  stress  of  a  square  plate  with 
similar  boundary  conditions  can  be  described  as: 
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The  critical  buckling  load  can  thus  be  normalized  with  respect  to  the 
tensile  yield  load  as: 


b. 


For  tha  mild  ataal  aquara  plataa  lnvaatlgatad  by  Ritohla  and  Rhodai, 
at turning  a  yiald  itrangth  of  30  kti,  tha  normalizad  critical  buckling 
load  la  roughly  0.106. 

Tha  analyaia  davalopad  by  Ritchia  and  Rhodai  datcribaa  a  normalizad 
ultimata  load  of  roughly  0,30,  with  raapact  to  tha  tanalla  yiald  load, 
for  thair  aquara  plataa.  Both  tha  critical  atraaa  and  tha  affactlva 
width  ara  atrongly  affaotad  by  ralatlva  plata  thleknaaa.  Tha  affactlva 
width  conoapt  datcribaa  a  normalizad  ultimata  load  for  an  unparforatad 
aquara  plata,  of  tha  a  am  a  dlmanaiona  and  matarlal  propartiaa  aa  taatad 
by  Ritohla  and  Rhodaa,  of  roughly  0.34  with  raapact  to  tha  tanalla  yiald 
load. 

Tha  thaoratioal  ralationthlp  batwaan  buckling  load  and  hola  aiza  in 
a imply  aupportad  plataa  loadad  by  both  uniform  traction  and  uniform 
dlaplacamant  la  aupportad  by  tha  axparlmantal  raaulta.  Tha  axparlmantal 
raaulta  of  othar  reaaarehara  for  uniform  traction  buckling  of  a  aquara 
plata  ara  plottad  in  Flgura  3.7  along  with  tha  analytic  aolution 
davalopad  by  Ritchia  and  Rhodaa.  Tha  axparlmantal  raaulta  for  uniform 
dlaplacamant  buckling  of  a  aquara  plata  ara  plottad  in  Flgura  3.8,  along 
with  tha  analytic  aolutiona.  Tha  changa  in  buckling  moda  in  ractangular 
plataa  from  two  to  thraa  half-wavaa,  whan  tha  hola  dlamatar  axcaada  0.33 
tlmaa  tha  plata  width,  waa  obaarvad  in  tha  axparlmantal  invastlgatlon. 
Tha  acattar  in  axparlmantal  raaulta,  according  to  Ritchia  and  Rhodaa, 
waa  dua  to  tha  larga  initial  imparfactlona  in  tha  plataa.  Tha  ultimata 
uniform  dlaplacamant  loada  for  both  aquara  and  ractangular  plataa  ara 
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presented  In  Figure  3.6  and  indicate  good  correlation  between  the 
analytic  solution  and  experimental  results.  Tha  ultimata  strength  of  a 
thin  plate  is  essentially  unaffected  by  the  presence  of  small  holes,  and 
is  roughly  equal  to  three  times  the  critical  buckling  load. 

With  regards  to  application,  however,  the  experimental  results 
Indicate  that  tha  elastic -plastic  ultimate  strength  analysis  developed 
by  Ritchie  and  Rhodes  is  slightly  conservative.  The  observed  normalized 
ultimate  load  for  a  thin  unperforated  square  plate  ranged  from  0,34  to 
0.35,  which  is  most  accurately  described  by  applying  the  effective  width 
concept.  It  is  suggested  that  the  Ritchie  and  Rhodes  analysis,  which 
accurately  predicts  the  loss  of  strength  associated  with  increased  hole 
size,  be  "shifted"  to  coincide  with  the  effective  width  solution  for  an 
unperforated  plate. 
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Figure  3.1:  Real  stress  distribution  and  sffsctivs  plats  width  dua  to 
in* plana  compressive  loading  [IS] 


Figure  3.2:  Graph  of  theoretical  uniform  traction  (constant  stress) 
and  uniform  displacement  (strain  buckling)  buckling  loads  versus  cutout 
ratio  for  square  plates  [16] 
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Figure  3.5:  Graph  of  affective  width  versus  maximum  stress  for  various 
cutout  ratios  [16] 
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Figure  3.6:  Graph  of  theoretical  analysis  and  experimental  results  for 
ultimate  load  versus  cutout  ratio  for  square  and  rectangular  plates  [16] 


Figure  3.7:  Graph  of  theoretical  analyst*  and  experimental  resulti  for 
uniform  traction  (constant  stress)  buckling  load  versus  cutout  ratio  for 
a  square  plate  [16] 


Figure  3.8:  Graph  of  theoretical  analysis  and  experimental  results  for 
uniform  displacement  (constant  strain)  buckling  load  versus  cutout  ratio 
for  a  square  plate  [16] 
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4.  DISHING  07  PLATES  UNDER  IMPULSIVE  LOADING 


4.1  Introduction 

Shipboard  structural  members  are  deformed  by  tha  over- 
pressurization  associated  with  an  internal  explosion.  The  effects  of 
holing  produced  by  the  initial  shock  wave  upon  residual  strength  of  a 
flat  plate  have  been  investigated  in  Chapters  2  and  3.  The  quasi- static 
overpressurization  which  follows  the  blast  produces  lateral  deflections 
in  the  bulkheads  and  deckplataa  of  compartments. 

The  deformation  of  thin  flat  plates  subjected  to  large  impulsive 
loading  has  been  studied  extensively  since  World  War  Two.  Thorough 
reviews  of  the  historical  development  of  empirical  and  analytical 
studies  have  been  presented  by  several  authors  [21,22,23].  A  number  of 
researchers  have  developed  good  analytic  techniques .  In  structural 
elements  such  as  beams  and  plates,  bending  resistance  increases  along 
with  membrane  forces  while  the  deflection  is  much  less  than  the 
thickness  of  the  element.  As  the  deflection  Increases,  however, 
stresses  due  to  stretching  of  the  middle  plane  dominate  over  bending 
stresses.  At  this  stage  the  behavior  of  a  beam  is  comparable  to  that 
of  a  plastic  string.  Perrons  and  Bhadra  [24]  developed  this  concept 
with  a  simple  lumped  mass  at  the  midpoint.  Recently,  Nurick  [23] 
developed  a  highly  accurate  membrane  mode  analysis  based  upon  the  lumped 
mass  string  concept.  Excellent  empirical  relationships  between 
deformation  and  Impulse  are  developed  from  extensive  testing  conducted 
by  Nurick  and  others. 


4.2  Pvnanlc  Efflfitl  &£  AQ  IntllMl  ALCblUt 


Roughly  70  percent  of  the  energy  released  by  a  warhead  detonation 
produces  blast  effects.  The  critical  loading  associated  with  an 
explosion  la  the  high  pressure  discontinuity  at  the  shock  front.  The 
shock  wave  produces  deflection,  necking,  rupture  and  holing  of  plating 
close  to  the  point  of  detonation  due  to  the  high  velocity  gradients 
across  the  plate  surface.  Rupture  begins  at  the  point  of  a  stiffened 
plate  structure  nearest  the  explosion  and  can  extend  to  the  plates' 
boundaries.  The  shock  loading  is  predominately  normal  to  all  structural 
members . 

After  passage  and  decay  of  the  shock  wave  from  an  Internal 
alrblaac,  confinement  of  the  explosion  gases  leads  to  a  pressure  loading 
which  has  a  long  duration  compared  to  the  response  time  of  typical 
stiffened  structures.  There  is  no  simple  relation  between  explosive 
mass  and  pressure  impulse  [  23  ] .  As  a  result  of  damage  to  bulkheads , 
decks,  and  overheads  surrounding  the  burst  space,  the  gas  overpressure 
is  vented  to  the  atmosphere  and  adjacent  spaoes  and  is  relieved  as  the 
confinement  volume  increases.  This  "quasi- static"  pressure  loading  is 
governed  by  the  confinement  volume.  Figure  4.1  shows  a  typical  trace  of 
the  pressure  versus  time  history  recorded  at  a  point  on  the  wall  of  a 
partially  closed  chamber  following  an  internal  explosion.  Pressure 
peaks  are  caused  by  repeated  reflections  of  the  blast  wave  off  the 
chamber  wall.  The  impulse  Imparted  by  the  pressure  peaks  is  small 
because  of  their  short  duration,  and  produces  negligible  damage  effects 
on  Internal  structure  and  equipment  [6,25]. 
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4.3  Mnahunt  Mod*  Flnlf  Element  An*lvli 

Membrane  mod*  approximations ,  In  which  bending  it  ignored,  have 
been  applied  to  impulsively  loaded  thin  plates  by  several  authors  who 
all  assumed  a  fixed  mode  shape.  Nurick  [23],  however,  has  developed  a 
modal  which  considers  lateral  displacements  in  addition  to  transverse 
deflections,  and  which  computes  the  mode  shape  at  each  time  step. 

The  string  concept  forms  the  basis  of  the  Nurick  model.  The 
string  is  sub-divided  into  several  sections  with  associated  lumped 
masses,  and  the  acceleration  is  assumed  to  be  a  linear  function  of 
the  velocity.  As  developed  by  Perrone  and  Bhadra  [24],  the 
acceleration  is  the  second  differential  of  displacement  resulting  in  a 
first  order  nonlinear  differential  eqt  tion  which  is  solved  by  numerical 
integration.  Transverse  and  lateral  deflections  are  incorporated 
into  the  equation  such  that  the  trajectory  of  every  point  on  the 
surface  remains  normal  to  the  plate  surface. 

The  material  is  assumed  to  behave  in  a  rigid  vlsco -plastic 
manner.  This  can  ba  described  by  the  empirical  relationship  developed 
by  Bodnar  and  Symonda  [26]: 


where  «  is  the  strain  rata,  D  and  n  are  material  constants.  Manjolne 
[27]  investigated  the  Influence  of  temperature  and  rate  of  strain  on  the 
yielding  of  mild  steel.  Bodner  end  Symonds  applied  Manjolne 's  work  to 
their  empirical  relationship  in  order  to  define  the  material  constants. 
For  mild  steel  the  material  constants  D  and  n  are  40  a'^  and  5 
respectively. 

The  deformed  geometry  of  the  string  can  be  defined  by  a  modeshape 
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4 ,  which  varies  with  apace  and  time.  The  transformed  displacement  is 
written  as: 

w  (r i t)  -  d  (r,t)  w  (t)  (4.2) 

y  f  0 

where  (r,t)  are  the  coordinates  of  space  and  time  and  w  (t)  is  the 

Q 

central  displacement.  It  is  further  assumed  that  the  transverse 
velocity  and  acceleration  fields  at  time  t  can  be  differentiated 
without  taking  into  account  the  variation  of  the  mode  shape  4  with  time, 

which  is  known  as  an  Instantaneous  or  stationary  mode  assumption. 

< 

The  string  concept  is  applied  to  plates  by  discretizing  the  plate 
into  contours  along  which  4  constant.  The  mode  shape  4  1*  assumed  to 
vary  linearly  with  distance  between  contour  levels.  Symmetry  dictates 
square  oontours  for  the  square  plates,  as  shown  in  Figure  4.2  .  This 
symmetric  contour  is  justified  by  the  deformation  shape,  shown  in  Figure 
4.3,  produced  by  explosive  impulse  loading. 

It  is  possible  to  write  (t  +  1)  dynamic  equations  for  the  lumped 
masses.  Nurick  developed  an  algorithm  to  solve  the  dynamic  equations. 
The  initial  central  velocity  is  computed  from  the  initial  total  impulse 
I,  which  is  assumed  to  impart  a  uniform  velocity  to  the  plate.  The 
initial  mode  velocity  is  then  obtained  from  a  generalised  momentum 
balance.  The  stiffness  of  ths  plate  is  assumed  to  be  sero  for  the  first 
time  step,  and  the  iterative  solution  commences  with  the  seoond  time 
step.  The  equations  are  highly  non* linear,  and  are  solved  using  a 
Newton- Raphson  technique.  The  Iterative  procedure  continues  until  a 
prespecified  tolerance  is  met  in  order  to  continue  to  the  next  time 
step,  and  the  forward  integration  is  terminated  when  w  equals  zsro. 
Convergeance  of  the  algorithm  suggests  that  the  final  shape  is  reached 
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within  100/is. 


4.4  Exairimantil  Rtiulti  And  Empirical  Antlyali 

The  experimental  investigation  of  Nurick  [23]  oomiatad  of  a  series 
of  tests  on  circular,  square  and  rectangular  plates  to  collect  a  range 
of  deflection  information  on  plates  subjected  to  high  impulsive  loading. 
Plate  specimens  were  cut  from  1.6mm  thick  cold  rolled  steel  and  were  not 
annealed.  Plate  dimensions  were  nominally  100mm  diameter  circular,  89 
by  89mm  square,  and  113  by  70mm  rectangular.  Metabel  was  used  aa  the 
sheet  explosive  for  the  experiments,  and  was  placed  on  a  16mm  thick 
polystyrene  pad  to  attenuate  the  shock  transmitted  to  the  plate,  provide 
a  uniform  impulse,  and  prevent  spallation  of  the  specimen.  The  specimen 
plate  was  clamped  between  two  holding  plates  attaohed  to  a  ballistic 
pendulum  used  to  measure  the  Impulse.  A  light  Interference  method 
was  used  to  measure  the  deflection* time  history.  The  final 

mid-point  deflection  was  measured  by  means  of  a  vernier  and  also  by  a 
reflex  metrograph. 

The  range  of  experimentation  performed  by  Nurick  was  impressive. 
Experimentation  by  previous  researchers  did  not  extend  beyond 
deflections  of  six  to  seven  plate  thicknesses  for  circular  plates  and 
eight  to  nine  thicknesses  for  quadrangular  plates.  Nurick  performed  a 
total  of  140  impulse  loading  experiments  with  deflections  ranging  up  to 
twelve  plate  thicknesses. 

The  analytic  prediction  method  developed  by  Nurick  compares 
favorably  with  impulse  loading  experimental  results  for  circular, 
square,  and  rectangular  plates.  The  range  of  analytic  predictions  and 
empirical  results  is  shown  in  Figure  4.4  for  circular  plates  and  Figure 
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4.5  for  square  plat**.  Although  most  of  th*  prediction  method*  can  do 
an  adequate  Job  predicting  deflection*  with  email  d*fl*ction-to- 
thicknes*  ratios,  they  display  tendencies  to  highly  over  or  under 
estimate  deflections  with  large  ratio*.  These  over  and  under  estimate* 
aro  based  upon  predictions  which  assume  a  fixed  mode  shape. 

For  circular  plates,  the  best  agreement  between  predictions  and 
experimental  results  was  achieved  by  Nurlck,  Llppman  [28],  Perrons  and 
Bhadra  [29],  and  Symonds  and  Wierrbicki  [22],  Th*  closed-form  solution 
for  circular  membranes,  developed  by  Symonds  and  Wlerzblckl,  is  a  simple 
but  accurate  approximation  method  which  disregards  bending.  For 
Impulses  which  blanket  the  plat*  surface,  their  solution  can  be 
simplified  as: 


w 


of 


.212  1 

gw* 


(4.3) 


where  w  is  the  final  deflection  at  th*  center  of  the  plate,  r  is  the 

0  a 

plat*  radius,  h  is  the  plat*  thickness,  and  p  is  th*  material  density. 
For  quadrangular  plates,  the  best  agreement  between  predictions  and 
experimental  results  was  achieved  by  Nurlck,  and  Baker  [30].  Using  an 
energy  method,  Baker  assumed  that  the  energy  Imparted  to  th*  structure 
could  be  equated  to  plastic  work  which  he  then  equated  to  th*  initial 
kinetic  energy  for  impulsively  loaded  structures.  Baker  assumed  that 
the  force  or  pressure  pulse  was  applied  in  a  time  short  enough  for  th* 
structure  response  to  be  only  a  function  of  the  Integral  of  th*  pulse, 
rather  than  th*  details  of  its  time  history. 

It  is  interesting  to  compare  and,  if  possible,  utilize  th* 
historical  data  available  from  different  impulse  loading 
experiments.  A  dimensionless  damage  number  was  defined  by  Johnson  [31] 
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where  v  the  impact  velocity  and  a  is  ths  damage  stress  (mean  flow 

O 

stress),  Table  4.1  shows  the  regime  of  damage  as  a  function  or  damage 
number,  and  shows  the  results  of  some  circular  plate  experiments.  The 
Johnson  damage  number  was  modified  by  Nurick  in  order  to 
non-dimenslonalize  dissimilar  geometries.  A  modified  damage  number  for 
circular  plates  can  be  written  as: 


9  - 

0 


*  r  ha  (fi  cr#) 


1/3 


<4. 5a) 


and  a  modified  damage  number  for  quadrangular  plates  can  be  written  as: 

.  _  I 


ha  (a  b  p  <r#)l/a 


(4. 3b) 


where  a  and  b  are  breadth  and  length  respectively  of  a  quadrangular 
plate,  and  o  is  given  the  value  of  the  static  yield  strength  ov. 

0  I 


Table  ft.l:  Typical  Johnson's  Number  Damage  Regime 
(mild  steel  at  room  temperature  impacted  by  a  projectile)  [31] 


ft/sec 

regime 

2.5 

10*9 

f  quasi* static 
\  elastic. 

25 

10** 

plastic  behaviour 

250 

10* 1 

slow  bullet  speeds 

2  500 

10  1 

extensive  plastic 
deformation 

25  000 

10  * 

hypervelocity  impact 

For  plates  of  equal  area  and  thickness  subjected  to  equal  total 
Impulses ,  the  midpoint  deflection  of  the  circular  plate  exceeds 
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that  of  the  square  plats,  which  In  turn  exceeds  that  of  tha 
rectangular  plate.  A  dimenaionleie  plot  of  deflection-to-thieknass 
ratio  veraua  modified  damage  number  8  for  quadrangular  plates  is 

9 

presented  in  Figure  4.6.  A  least  squares  correlation  analysis  of 
the  experimental  results  indicates  that  Murick's  experimental  results 
are  consistent  with  previous  data,  and  that  a  relationship  exists 
between  the  impulse  a  given  plate  is  exposed  to  and  the  deformation  it 
experiences . 

Based  upon  the  extensive  range  of  consistent  experimental 
results,  it  is  possible  to  develop  an  empirical  solution  to  the 
deflection  problem.  A  least  squares  analysis  of  the  olrcular  plate 
experimental  data  yields: 

-  0.423  +  0.277  (4.6a) 

'  hare  the  number  of  data  points  was  109  with  a  correlation  coefficient 
of  0.974.  A  least  squares  analysis  of  the  quadrangular  plate 
experimental  data  yields: 

[■Tv1]*  -  0.236  +  0.001  (4.6b) 

where  the  number  of  data  points  was  156  with  a  correlation  coefficient 
of  0.984.  The  empirical  relationships  can  be  simplified  to: 


[H 

-  0.425  8 

e 

(4.7a) 

r  w  _  -i 

[■*"], 

-  0.236  8 

q 

(4.7b) 

as  a  useful  "engineering"  approximation. 

The  excellent  correlation  of  the  empirical  formulas  Indicate  that 
they  can  be  readily  applied  to  engineering  analyses.  The  linear  nature 
of  the  closed- form  deformation  solution  developed  by  Symonds  and 
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Wlarzbickl,  and  othara,  supports  tha  uaa  of  a  linaar  empirical 
ralatlonahip  basad  upon  a  modlfiad  Johnson's  damaga  number.  Extanslva 
deformations  can  produca  taaring  along  tha  adgas  of  tha  plata  or  holing 
In  tha  center  of  tha  plata.  Tha  maximum  non- taaring 

deflation- to -thlokneas  ratios  during  Nuriek'a  tasting  wars  11.60  and 
11.97  for  circular  and  quadrangular  plates  respectively.  This  taaring 
along  tha  clamped  boundaries  indicates  that  both  analytical  and 
empirical  prediction  methods  are  limited  to  modlfiad  damaga  values  of  9 
&  25  for  circular  plates  and  6 at  50  for  quadrangular  plates. 
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Figure  4.3  Final  mod*  profil*  and  contour  plot  of  *  typical  explosive 
tast  for  a  squat*  plat*  [23} 


IMPULSE  (Pi) 


Figure  4.4:  Graph  of  analytic  solutions  and  experimsntal  results  for 
oeflection- to- thickness  ratio  versus  impulse  for  circular  plates  [23] 


5.  STRENGTH  07  A  DISHED  PLATE 


5.1  Introduction 

Despite  the  large  lateral  deformations  suffered  by  shipboard 
bulkheads  and  deckplates  as  a  result  of  impulsive  loading,  the  plating 
can  continue  to  support  normal  service  loads.  The  dishing  of  thin  flat 
plates  into  a  shell  was  related  to  such  parameters  as  material 
properties  and  explosive  Impulse  in  Chapter  4.  By  examining  large  scale 
elastic -plastic  deformation  mechanisms  it  is  possible  to  describe  the 
buckling  behavior  of  a  thin  shell. 

An  out-of-flatness  plate  with  a  small  initial  deflection  does  not 
react  to  in-plane  loading  like  an  ideal  flat  plate.  The  effect  of 
initial  out-of-plane  deformations  is  generally  to  reduce  the  peak  load 
of  the  plate  and  bias  the  mode  of  failure  towards  the  initial  shell 
shape.  The  effect  of  initial  imperfections  upon  the  ultimate  strength 
of  plates  offers  some  insight  into  the  effect  of  large  scale  dishing. 
The  residual  strength  of  a  thin  shell  with  large  initial  deformation  (on 
the  order  of  5  to  15  plate  thicknesses)  is  first  analyzed  as  a  simple 
extension  of  the  concept  of  initial  imperfections.  An  extenslonal 
collapse  mode  analysis  is  then  developed  as  a  simple,  closed- form 
solution  which  accounts  for  actual  deformation  mechanisms. 

5.2  Efrlilfill  BmIi  fi£  Deformation  Mechanisms 

The  yield  load  of  a  plate  and  the  elastic-plastic  buckling  load 
provide  two  limits  to  plate  strength.  The  onset  of  plasticity  is 
affected  by  both  the  in-plane  and  plate  bending  stress  components. 


Plasticity  becomes  more  significant  as  the  initial  out-of-plane 
imperfections  become  larger  causing  surface  and  membrane  yield  to  occur 
earlier.  The  Interaction  of  membrane  and  bending  atresses,  with  initial 
imperfections,  produces  an  early  loss  of  stiffness  and  reduction  in  peak 
load  balow  the  Ideally  planar  elastic  buckling  load  [17], 

Although  the  ultimate  strength  of  a  plate  with  very  slight  initial 
imperfections  remains  raasonably  close  to  that  of  an  ideal  flat  plate, 
the  effects  of  larger  imperfections  (on  the  order  of  the  plate 
thickness)  can  be  significant  [17],  Von  Kdrmdn  daveloped  the  effective 
width  concept  (presented  in  Chapter  3)  to  describe  the  stress 
distribution  produced  in  a  simply  supported  plate  subjected  to  uniaxial 
compreaaive  loading.  The  Influence  of  geometrical  imperfections  has 
been  incorporated  into  a  von  Kdrmdn  analysis  by  several  authors ,  as 
shown  in  Table  5.1.  Initial  deflections  of  the  plate  lead  to  a  marked 
reduction  of  the  peak  load,  especially  for  thin  plates  with  reference 
slenderness  Ap  less  than  1.5.  From  Eq.  3.2,  the  reference  slenderness 
is; 


X  ■ 

P 


b  1,05 

K  /TT 


(5.1) 


Rectangular  plates  restrict  the  failure  mode  and  generally  have  greater 
stiffness  than  square  plates,  but  suffer  a  steeper  post-buckling 
behavior  (17). 

In  order  to  understand  the  true  collapse  behavior  of  a  shell 
subjected  to  uniaxial  compressive  loading,  the  effect  of  out-of-plane 
deformations  on  the  in-plane  stress  state  of  the  plate  must  be 
considered.  In  the  large  deflection  regime,  in  general,  in-plane  and 
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out -of -plan*  displacements  end  hence  stresses  end  strelns  interect.  As 
e  result  of  dishing,  the  length  of  the  plete  center  line  becomes  greeter 
then  the  boundery  length  end  membrene  tensions  develop  in  both  plete 
directions.  Membrene  ection  becomes  more  end  more  dominent  es 
out-of-plane  deflections  become  larger  [17]. 

Table  5.1:  Empirical  relationships  to  describe  the  effect  of 
smell  initial  imperfections  upon  effective  plete  width  [17] 


Winter , 

bt 

r  “ 

M1- 

for  Ap  fc  0.67 

Faulkner , 

bt 

E~  ' 

1.05  f. 

it  L1 

for  Ap  fc  0.55 

Gerard, 

bt 

0.82 

ET  “ 

.  0.19 

A* 

5.3  km  Initial  Imgis, flatten  Analxite 

Extensive  research  has  been  conducted  in  the  lest  twenty- five  years 
on  the  ultimate  load  behavior  of  box  girder  structures  [32,33,34],  This 
research  has  followed  a  series  of  structural  failures  involving  merchant 
ships  (large  ore  carriers  in  rough  seas  near  Japan)  [33]  and  large  box 
girder  bridges  in  Europe  [34].  Frieze,  Dowling  and  Hobbs  [32]  have 
developed  a  large  deflection  elastic-plastic  analysis  to  examine  the 
behavior  of  plates  in  compression. 

The  analytic  method  developed  by  Frieze  et  al.  uses  a  modified  von 
KdrmAn  analysis  to  account  for  material  non-linearity  and  the  effects  of 
residual  stresses  end  initial  out-of-plane  deflections.  The  governing 
equations  were  solved  by  computer  processing  using  an  iterative  finite 
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difference  procedure.  Frieze  et  el.  performed  e  parametric  study 
examining  the  collapae  behavior  of  a  unlaxially  loaded  plate.  This 
analyaia  indicate!  that  the  influence  of  expect  ratio  (a/b)  ia  dependent 
upon  both  the  magnitude  of  the  initial  deformation  and  the  slenderness 
ratio.  Although  the  in-plane  restraint  along  unloaded  edges 
(constrained  or  unrestrained)  has  little  effect  on  the  behavior  of 
stocky  plates,  the  analysis  indicates  that  slender  plates  with 
constrained  edges  retain  their  strength  better  than  with  unrestrained 
edges.  Figure  5.1  shows  that  this  is  especially  true  for  plates  with 
large  initial  deformations.  Slender  plates,  in  particular,  suffer  from 
the  weakening  effects  of  initial  deformations  which  generally  dominate 
the  effects  of  residual  stresses. 

/ 

The  residual  strength  of  a  thin  shell  can  be  extrepolated  from 
theoretical  and  empirical  analyses  which  have  been  performed  for  modest 
initial  deformations.  It  is  not  possible  to  directly  apply  the  analyses 
because  shipboard  structural  plating  is  very  thin.  The  analyses  can  be 
extrapolated  by  applying  the  von  KArmAn  effective  width  concept  to 
describe  the  loss  of  strength  for  unperforatod  flat  plates  with 
increasing  slenderness.  The  theoretical  analysis  of  Frieze  et  al.  is 
summarized  by  Figure  5.2,  which  presents  six  different  curves  of 
ultimate  load  versus  slenderness  fi  for  different  levels  of  initial 
dishing.  The  slenderness  value  fi  is  defined  as: 

'-Hr*  <5,2> 

which  is  roughly  twice  the  plate  reference  slenderness  Ip  for  a  square 
plate.  The  slenderness  ratio  for  shipboard  structural  plating, 
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disregarding  tha  affect  of  stiffeners,  ranges  roughly  from  9  to  IS. 
Ueda  at  al.  [35]  also  examined  the  effects  of  initial  deformations  and 
residual  stresses  on  the  compressive  strength  of  square  plates,  and 
conducted  a  modest  program  of  experimentation.  Large -deflection, 
elastic-plastic  analyses  were  performed  using  an  incremental  finite 
element  technique.  Usds' s  theoretical  and  experimental  research 
supports  Freize's  analysis,  as  shown  in  Figure  5.3,  and  by 
extrapolation  suggests  that  ths  normalized  ultimate  load  would  decrease 
for  a  thin  shell  with  large  initial  deflections  by  roughly  10  percent 
from  that  of  a  flat  plate  with  similar  slenderness.  The  plate  critical 
buckling  curve  is  included  in  both  Figures  5.2  and  5.3  for  comparison 
purposes.  Even  with  large  initial  imperfections,  the  ultimate  strength 
of  slender  plates  is  greater  than  the  critical  buckling  load. 

5.4  Extenslonel  Collapse  Mode  Analysis 

A  structure  can  absorb  crushing  work  through  local  plastic  folding 
with  bending  and  extans lonal  deformation.  Minorsky  [36]  examined  the 
ship  collision  problem,  with  an  analysis  of  smpirlcal  data  from  several 
dozen  reported  casualties,  and  found  that  plastic  work  dissipation  is 
linearly  related  to  the  volume  of  damaged  members.  Although  useful  for 
studying  major  collisions,  Minorsky's  method  is  not  valid  for  structural 
failure  of  a  local  nature.  The  theory  of  extenslonel  collapse  of  plate 
elements  was  developed  by  Hayduk  and  Wierzblcki  [37]  and  others  in  the 
early  1980 's  to  examine  the  crush  behavior  of  Individual  members  during 
crushing.  Extenslonel  collapse  mode  analysis  is  described  here  in  order 
to  determine  the  average  load-carrying  capability  in  structural  plating 
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crushed  beyond  tha  ultimata  load. 

Ships  ara  generally  oonstruotad  of  plate-like  structures,  aa  shown 
in  Figura  5.4.  During  dafornation,  tha  adgas  of  tha  structure  can  ba 
considarad  hinga  lines,  as  can  tha  horizontal  limits  of  tha  buckling 
region.  Crushing  work  is  dissipated  through  plastic  banding  of  adjacent 
plates  about  the  horizontal  hinge  lines  connecting  them  and  through 
shape  distortion  of  each  plate.  The  energy 'balance  postulated  by 
Wierzbioki  and  Abramowiez  [38]  indicates  that  zones  of  bending  and 
extansional  deformation,  which  are  restricted  to  a  small  fraction  of  tha 
total  area  of  the  structure,  dissipate  essentially  all  of  tha  work  dona 
by  the  applied  load  during  tha  crushing  deformation  mode. 

In  the  propagating  hinge  line  deformation  mode  (Mode  X),  extension 
is  limited  to  the  neighborhood  of  the  hinge  line.  Displacement  fields 
must  be  described  by  a  continuous  and  smooth  funotiqn  across  a 
propagating  hinge  line.  The  quas 1 - lnaxtens ional  eollapse  mode  of  a  L- 
shaped  angle  element  is  shown  in  Figure  5.5.  Wierzbioki  and  Abramowiez 
described  the  individual  dissipation  energies  associated  with  the 
deformation  of  a  toroidal  shell,  horizontal  hinge  line,  and  Inclined 
hinga  line  through  complete  folding  of  the  section.  In  the  basic 
folding  mechanism  of  a  L- shaped  angle  element  there  are  eight  horizontal 
hinge  lines,  two  inclined  hinge  Unas,  and  one  toroidal  shell  section, 
The  resulting  theoretical  crushing  strength  of  a  L-shapad  element  is 
thus  given  by: 

r  * 12-04  J[  is )  <5,3) 

where  is  tha  mean  crushing  load  and  is  the  fully  plastic  yield 
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moment  defined  u : 


Short-angle  columns  commonly  fold  up  and  fall  in  an  aaymmatrlc  mods 
uniats  vary  largo  imparfoctions  art  lntroducad  in  Cha  shapa  of  tha 
symmatric  modas  [37]. 


In  tha  stationary  hinga  lina  daformation  moda  (Mode  II)  tha  adgas 
of  tha  structura  remain  stationary  whlla  tha  slda  plataa  undargo 
rigid-plastic  distortion.  Plata  distortions  consist  of  althar  axtansiva 
in-plane  oomprassion  or  in-plana  extension.  Hayduk  and  Wlarcblckl 
assumad  that  thera  is  a  linaar  variation  of  tha  valooity  flald  In  tha 
horlsontal  y-dlraotlon,  as  shown  In  Figura  5.6,  and  davalopad  tha 
individual  dissipation  anarglas  assooiatad  with  tha  daformation  of  a 
plata,  stationary  horizontal  hinga  and  stationary  incllnad  hinga  through 
complata  folding  of  tha  saotlon. 

Hayduk  and  Viarzbieki  appliad  tha  axtanslonal  collapsa  moda 
analysis  to  a  x-shapad  cruclfom  structural  alamant.  A  cruciform  Is 
constructed  of  four  flanges  which  are  rigidly  joined  along  tha  initially 
vartloal  adga.  A  cruciform  alamant  reasonably  simulates  tha  boundary 
conditions  on  a  plata  in  a  typical  shipboard  compartment.  Within  tha 
folding  zona,  tha  joint  line  becomes  an  incllnad  hinga  ,  fixed  in  tha 
material,  In  Mode  II  crushing  of  a  cruciform  there  are  eight 
trapezoidal  distortion  sections,  sixteen  stationary  horizontal  hinge 
Unas,  and  four  stationary  incllnad  hinge  Unas.  An  upper  bound  to  tha 
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crushing  strength  of  a  cruciform  alamsnt  can  thus  ba  described  as: 


5T  "  20l°5 
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(5.6) 


This  least  upper  bound  is  based  upon  a  minimum  energy  analysis  which 
indicates  that  folding  will  occur  with  a  zone  of  half -wavelength  H: 


(5.7) 


The  upper  bound  analysis  also  indicates  that  the  energy  dissipated  at 
the  inclined  hinge  lines  is  negligible  for  a  sufficiently  large 
length* to- thickness  ratio. 

A  lowor  bound  on  the  crushing  strength  of  a  cruciform  can  be 
similarly  derived  by  noting  that  a  continuous  Joint  along  the  vertioal 
edge  of  a  cruciform  increases  the  column  strength  as  oompared  to  a 
column  composed  of  two  unjoined  L* shaped  elements.  From  the  Mode  I 
analysis  of  a  single  L* shaped  element,  a  lower  bound  to  the  crushing 
strength  of  a  cruciform  can  ba  described  as: 


ft  *  24‘08  /[i] 


(5.8) 


Careful  inspections  of  several  crushed  cruclforms,  such  as  the  one  shown 
in  Figure  S.7,  reveal  that  the  actual  collapse  mechanism  is  a 
combination  of  the  two  basic  deformation  modes.  Hayduk  and  Wiarzblcki 
developed  a  modified  Mode  II  model,  and  described  the  mixed-mode 
crushing  strength  of  a  cruciform  as: 
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(5.9) 


The  Hayduk  and  Wiarzblcki  extenslonal  collapse  mode  solutions  for 
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cruciform  failure  are  baaed  upon  complete  folding  of  the  atructure, 

Thia  extenaional  collapae  mode  analysis  of  a  cruciform  can  be 
applied  to  the  reaidual  atrength  of  a  diahed  plate.  The  typical 
load-deflection  characteriatica  of  a  cruahed  cruciform,  ahown  in  Figure 
5.8,  indicate  that  the  ultimate  strength  of  a  cruciform  with  Initially 
flat  aide  platea  la  rapidly  reached  during  uniaxial  loading.  The 
post-buckling  and  folding  behavior  after  the  ultimate  atrength  la 
reached  "reinforces11  ths  cruciform  and  allows  it  to  continue  to  support 


loading  during  repeated  folding.  Since  the  collapae  mode  analyla 
indicates  that  folding  will  occur  along  a  zone  of  width  2H  defined  by 
Eq.  5.5,  dishing  may  not  significantly  reduce  a  plate's  ultimata 
atrength  because  the  diahed  shape  does  not  encourage  the  folding. 
Although  dishing  may  reduce  a  plate's  ultimate  atrength,  it  la 
postulated  that  the  mean  crushing  load  Which  the  structure  can 
support  is  little  affected  by  the  dishing. 

The  extenaional  collapse  mode  analysis  equations  can  be  rewritten 
to  provide  normalized  loads  for  shell  crushing.  The  mean  crushing  load 
la  normalized  with  respect  to  the  tensile  yield  load,  By  assuming  that 
the  strength  of  a  single  dished  plate  can  be  considered  equal  to 


one-half  that  of  a  cruciform,  the  upper  bound  can  be  expressed  as: 


a 

e 


2.51  h  [b/2h)1/l 
b 


(5.10a) 


the  lower  bound  can  be  expressed  as: 


o  3.01  h  [b/2h] 


(5.10b) 


and  the  mixed-mode  solution  can  be  expressed  as: 
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?m  ao  2.87  h  {b/2h)l,i  +  0.36  h  (b/2h) 
* 


2/9 


(3.10c) 


*  b 

A  modern  warship's  structural  plating,  disregarding  the  effect  of 
•tlf fenara ,  hea  a  typical  length- to- thickness  ratio  (b/h)  on  the  order 
of  200  to  300.  From  the  upper  end  lower  bound  analysis  for  a  warship's 
structural  placing,  the  normalised  mean  crushing  load  ranges  from  0.05 
to  0.13. 

The  axtensional  collapse  mode  analysis  can  be  applied  to  a  holed 
shell.  Since  the  width  of  the  local  folding  sons  is  much  smaller  than 
the  height  of  the  typical  cruciform,  as  shown  in  Figure  5.9,  the  "wing 
areas"  can  be  considered  ineffective  in  terms  of  both  axtensional  and 
bending  energy.  By  assuming  that  the  fold  is  made  through  the  hole 
region,  the  affective  plate  half -width  can  be  reduced  by  the  radius  of 
holing,  For  Eqs.  5.10a  •  5.10c  only  the  half-width- to- thickness  ratio 
(b/2h)  in  the  numerator  is  redueed  by  holing.  The  axtensional  oollapse 
mode  analysis  can  then  be  applied  using  the  reduced  half -widths. 


5 . 5  Comparison  With  Experimental  Results 

The  experimental  investigation  of  Hayduk  and  Wierzbickl  [37] 
consisted  of  a  series  of  tests  on  thin  copper  cruciform  specimens. 
Specimens  ware  constructed  of  two  L- angle  elements  of  0.014  inch  thick 
copper  sheeting  butted  along  the  bend  line  and  soldered  with  common 
60-40  lead-tin  solder.  Specimen  length  varied  from  2.0  to  4,0  in,  while 
plate  half-widths  varied  from  1.0  to  2.0  in.,  for  a  length -to -thickness 
ratio  which  varied  from  143  to  286.  The  unloaded  vertical  edges  of  the 
cruciform  were  clamped  at  the  top  and  bottom  to  prevent  specimen 
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twitting  during  loading.  During  quaii-itatic  loading  eha  tptcimtnt  wtrt 
routinely  cruthtd  from  40  to  60  ptrctnt  of  their  original  length. 

The  extent  ional  cruthing  mode  anti.;  rt  it  it  tupported  by  the 
experimental  retultt.  Figure  5.10  presentt  the  experimental  data  and 
the  upper,  lower  and  mixed-mode  tolutiont  for  oruthing  of  copper 
cruciformi  and  other  elementt,  The  experimental  data  for  cruciformt  are 
thown  by  open  circlet.  The  tcatter  of  experimental  data  lndicatei  that 
the  cruciformt  failed  according  to  a  varying  combination  of  collapae 
mechanitmt.  The  collapae  mechanitmt  are  tentitive  to  imperfectiont. 
Indication!  of  fracture  at  the  joint  in  tome  tett  piecet  explain!  tome 
of  the  tcatter  below  the  mixed-mode  eolution.  Deipite  the  tcatter  of 
data,  all  experimental  retultt  for  cruciformt  lie  within  the  upper  and 
lower  bound  tolutiont  to  normalised  mean  cruthing  load.  Amdahl  [39] 
performed  a  timllar  teriet  of  tettt  on  larger  iteel  cruciformt  whioh 
alto  tupport  the  validity  and  utefulneat  of  extemlonal  cuthlng  mode 
analytli, 
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SLENDERNESS  C/?3 

Figura  5.3:  Graph  of  thaoratleal  analyala  and  axparimantal  raaulta  for 
ultimata  load  varsus  slandarnaaa  for  aquara  platai  including  tha  affacta 
of  Initial  daflactlona  and  voiding  raaidual  itraaaaa  [33] 
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Figure  5.5:  Basic  folding  mode  of  a  L- shaped  angle  element  with  a 
continuous  and  smooth  displacement  field  [37] 


Figure  5.6:  Plastic  distortion  of  a  rectangle  into  a  trapezoid  [37] 
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Figure  5.9:  Cruciform  with  holed  half -plates 
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Figure  5.10:  Graph  of  theoretical  analysis  and  experimental  results  for 
mean  crushing  strength  versus  half -plate  length- to- thickness  ratio  for 
cruciforms  [37] 
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Figure  5.9:  Cruciform  with  holed  half -plate* 
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Figure  5,10:  Graph  of  theoretical  analyaia  and  experimental  result*  for 
mean  crushing  strength  versus  half-plate  length- to -thickness  ratio  for 
cruciform*  [37] 


6.  MIXED  MODE  CRACK  PROPAGATION 


6.1  Introduction 

The  holt  edge  produced  during  the  fracture  and  taarlng  of  natal 
plataa  aasoolatad  with  an  lntarnal  alrblast  la  not  smooth.  Although  tha 
hola  can  reasonably  ba  consldarad  to  ba  circular  for  soma  forms  of 
analysis,  It  is  Important  to  examine  tha  affaot  of  a  raggad  and  torn 
adge  upon  rasldual  strangth.  Initial  cracks  of  various  sizes  will 
amanata  from  around  tha  rough  adga  of  tha  hola.  Crack  growth  undar 
loading  (monotonic  or  cyclic)  may  limit  tha  rasldual  strangth  of  tha 
structura. 

Tha  initiation  of  cracks  during  holing  and  ovarpraasurlsatlon 
assoclatad  with  tha  airblast  is  a  complax  dynamic  situation.  Tha 
daformad  plata  can  ba  oonsldarad  a  thin  shall.  Crack  initiation  and 
propagation  from  tha  torn  adgas  assoclatad  with  holing  and  fragmantation 
occurs  in  a  fully  plastic  or  elastic-plastic  deformation  field. 

Crack  propagation  ooours  as  a  mixture  of  three  limiting  modes,  as 
shown  in  Figure  6.1.  Mode  I  crack  propagation  Is  assoclatad  with 
in-plane  tension,  Mods  II  is  assoclatad  with  ln-plane  shear,  and  Mode 
III  is  associated  with  out-of-plane  shear.  Actual  crack  propagation  is 
rarely  so  distinct  (a.g.  slant  fracture  consisting  of  Modes  I  and  III 
can  ba  caused  simply  by  Mods  I  loading).  As  one  could  imagine  tearing  a 
page  out  of  a  telephone  book,  tha  metal  plating  will  ba  opened  by 
ln-plane  tension  and  torn  by  out-of-plane  shear.  Crack  propagation  from 
the  rough  edge  of  a  hola  in  a  thin  shell  will  be  caused  primarily  by 


6.2  Physical  Reel*  of  Pla«tle  Pratt tura 


A  fully  plastic  situation  consists  of  piano  strain,  plana  atraaa, 
or  intamadiata  nodas  which  occur  whan  tha  plastic  zona  is  largo 
ralativa  to  tha  plata  thicknass.  Tha  plastic  zona  will,  for  plana 
stress,  become  larga  and  spraad  through  tha  whola  crackad  section  if  the 
net  section  stress  is  of  tha  ordar  of  tha  yiald  strangth.  In  practica, 
linaar  alastie  fracture  mechanics  (LEFK)  does  not  apply  if  tha  fracture 
stress  is  greater  than  two-thirds  of  yiald  [40],  Plana  stress  requires 
tha  vanishing  of  one  of  tha  principal  components  of  stress,  and  tha 
vanishing  of  stress  gradients  in  tha  s  direction.  Plana  stress  occurs 
whan  tha  out -of -plana  dimension  is  far  thinner  than  that  required  for 
plane  strain,  as  in  a  thin  plata  or  shall. 

Tha  tearing  of  thin  sheets  of  ductile  metals  is  a  process  of 
oomblned  flow  end  fraotura.  Plasticity  precedes  and  accompanies  crack 
propagation.  Necking  may  occur  whan  the  plastic  zona  size  becomes 
larger  than  the  plate  thickness.  Stress  end  strain  distributions  in  the 
region  of  tha  tip  of  a  sharp  oraok  are  complex.  Stress  ohenges  are  vary 
large  within  one  plate  thickness  of  the  crack  tip,  and  a  state  of  plane 
strain  may  be  approached.  Plana  stress  behavior  will  be  exhibited 
beyond  roughly  four  plata  thicknesses  of  the  craok  tip  [10] . 

Crack  propagation  can  be  regarded  as  a  process  of  hole  nucleatlon, 
coalescence,  growth,  and  localization.  Localization  may  occur  first  in 
plana  stress.  Tipper  [41]  examined  sections  near  tha  fracture  surface 
cut  from  different  mild  steal  specimens.  In  ductile  steels,  hole 
nucleatlon  almost  certainly  occurs  by  the  drawing  away  of  the  metal  from 
non-metalllc  inclusions.  Separation  of  the  metal  around  inclusions 
occurs  at  some  distance  (1  or  2  hole  spacings)  ahead  of  the  crack  tip, 
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and  may  occur  with  stresses  aa  low  aa  yield.  The  holea  become  linked  by 
elongation  and  eventually  produce  macrofracture. 

6.3  Analytic  APPIttlChf  Lfi  Pintle  Fracture 

Plaatic  fracture  ia  a  complex  proceaa  whose  mechanisms  are  not 
fully  understood,  These  mechanisma  Include  hole  nuclaation, 
coalescence,  growth,  and  localisation  of  plaatic  flow  into  banda. 
Linear  and  non-linear  elastic  fracture  mechanica  are  relatively  well 
understood  for  crack  initiation  in  plane  strain,  but  do  not  apply  to  the 
loading  of  shells  and  subsequent  crack  growth.  Despite  considerable 
research  within  the  last  tan  years,  according  to  Broek  [40],  no 
generally  accepted  approach  to  elastic-plastic  fracture  mechanics  (EPFM) 
has  been  developed.  Kardomataas  and  KoCllntook  [42]  indicate  that 
ductile  fracture  occurs  by  hole  growth,  which  depends  on  the  combined 
history  of  strain,  stress,  and  rotation  near  the  orack  tip.  Fully 
plastic  fracture  mechanics  (FPFM)  can  thus  not  be  described  by  a  single 
measure  of  failure. 

Several  promising  approaches  to  EPFM  are  being  developed,  however, 
which  offer  some  insights  into  the  process.  These  approaches  generally 
consist  of  either  a  strass  field  analysis  or  an  approximate  global 
energy  analysis.  A  strass  field  approach  examines  the  stress  field 
around  the  crack  tip,  as  the  K  factor  does  in  LEFM.  An  energy  approach 
examines  the  energy  release  associated  with  crack  initiation.  Both 
approaches,  unfortunately,  have  their  limitations. 

Although  it  has  not  been  accepted  as  a  universal  fracture 
criterion,  the  J- integral  provides  a  means  to  determine  an  energy 
release  rate  for  both  LEFM  and  cases  where  plasticity  effects  are  not 
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negligible  [40].  The  J - integral  ie  defined  along  a  contour  around  the 
crack  tip  aa  the  change  in  potential  energy  for  a  virtual  crack 
extanalon.  The  U.S.  Navy  has  adopted  a  J-based  approach  to  fracture 
toughneaa  characterizationa  of  ductile  materlala  [43].  In  general, 
though,  the  J- integral  ia  unable  to  handle  large  amounts  of  stable  crack 
growth  accompanied  by  full  plasticity  because  unloading  occurs  along  a 
different  stress-strain  curve  than  during  loading,  and  also  because  of 
differing  degrees  of  triaxiality  with  different  geometries  and  low 
hardening.  McCllntock  [10]  further  emphasizes  that  an  energy  approach 
to  plastic  crack  propagation  is  global  in  nature  and  can  not  account  for 
the  different  fracture  process  zone  stress-strain  distributions  and  the 
local  mechanisms  of  crack  growth.  Craok  initiation,  however,  is  an 
especially  local  phenomenon  and  can  be  adequately  described  by  an  energy 
approach. 

Atkins  [44]  has  developed  a  rigid-plastic  fracture  analysis  which 
reasonably  accounts  for  and  simplifies  the  complex  nature  of  fully 
plastic  fracture  (FPFM)  in  certain  situations.  In  this  approach,  the 
relative  magnitude  of  flow  and  fracture  components  is  an  Important 
factor  In  determining  the  transition  from  plastic  flow  proceeding 
fracture  to  crack  initiation  and  propagation.  Although  this  approach  is 
limited  by  its  global  nature,  it  offers  a  useful  engineering  analysis  of 
FPFM,  and  is  presented  because  of  its  promise. 

A  limit  load  analysis  offers  a  new  approach  to  the  problem  of  fully 
plastic  crack  propagation.  By  selecting  realistic  displacement  fields 
it  may  be  possible  to  define  what  loads  cause  deformation  which  could 
tend  to  promote  crack  propagation. 
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6.4  Rigid. Plastic  Fracture  Anal Vila 


A  rigid-plastic  simplification  of  EPFM  mav  bo  applied  to  'a  variety 
of  sheet  tearing  problems.  Plastic  fracture  Involves  remote  flow,  which 
Atkins  [44]  incorporated  into  an  expanded  quaai.static  elastic  fracture 
energy  equation  in  terms  of  the  elastic  strain  energy  A,  the  fracture 
toughness  or  specific  (essential)  work  of  fracture  R,  the  total  plastic 
work  Wp,  the  residual  strain  energy  Ag ,  and  kinetic  energy  H  as: 

Pdu-dA  +  RdA  +  dWp  +  dA^  +  dH  (6.1) 

where  P  is  load,  u  is  load.point  displacement,  and  A  is  the  crack  area. 
In  the  specific  case  of  dynamic  tearing  of  ductile  metal  sheets  the 
limiting  crack  speeds  consistent  with  a  quasi-static  approximation  are 
usually  some  1  to  10  m/s  [44] .  The  kinetic  energy  term  may  be  equated 
to  zero  for  most  materials  and  cracked  geometries.  Cracks,  whether 
surrounded  by  elastic  or  plastic  fields,  tend  to  propagate  with  their 
own  crack  opening  angles.  The  rigid-plastic  fracture  equation  can  be 
further  simplified  in  the  case  of  extensive  plastic  flow,  where  it  may 
be  permissible  to  neglect  dA,  as: 

P  du  -  R  dA  +  dWp  .  (6.2) 
This  equation  has  found  success  in  Cotterell  and  Hal's  [45]  studies  of 
fracture  in  ductile  sheets,  and  in  Atkin's  studies  of  metal  forming 
fracture  such  as  cropping  and  can  opening.  These  studies  demonstrated 
that  geometry* independent  R  values  do  characterise  the  fracture  behavior 
and  that  plastic  flow  and  fracture  terms  may  be  separated  in  the 
rigid-plastic  fracture  aquation.  Atkins  cautions,  however,  that 
fracture  toughness  may  be  path  dependent  according  to  a  McCllntock-type 
analysis . 
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Atkins  and  others  have  examined  crack  propagation  due  to  various 
loading  modes.  Experiments  by  Mai  and  Cotterell  [46]  in  which  the 
banding  and  unbending  plastic  work  could  be  altered  Independently  of  the 
fracture  work  gave  roughly  the  same  values  for  fracture  toughness.  The 
experiment  used  two  metals  (low  alloy  steel  Lyten  and  a  5251  allumlnum 
alloy)  with  different  work  hardening  indices  in  both  two  and  three  1*1 
"trouser”  tests  to  alter  the  plastic  work  independently  of  the  fracture 
work.  This  demonstrates  that  the  flow  and  fracture  work  components  can 
sometimes  be  uncoupled  and  that  the  toughness  so  determined  is  a 
material  property,  and  indicates  that  the  rigid-plastic  equation  is 
valid  for  certain  applications. 

From  a  global  energy  approach,  crack  propagation  will  occur  if  the 
oombined  flow  and  fracture  associated  with  propagation  requires  a 
smaller  work  increment  than  flow  alone.  An  increment  of  plastic  flow 
may  be  written: 

dWp  -  d(r  V)  -  T  dV  +  V  dT  (6.3) 

where  i"  is  the  total  plastic  work  per  unit  volume  and  V  is  the  volume 
undergoing  plastic  deformation.  The  plastic  flow  aquation  governs  the 
external  work  required  prior  to  fracture  and  applies  with  Incremental 
changes  dV  and  dWp  since  there  is  new  material  entering  the  fracture 
process  zone  as  crack  propagation  occurs.  By  considering 
Increments  of  external  work,  and  the  additional  work  required  to 
propagate  cracks,  Atkins  was  able  to  describe  stable  crack 
propagation  as: 


where  n  is  the  number  of  propagating  splits  [44].  The  stability  of 
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crack  propagation  la  dapandant  upon  boundary  conditions  and  and  loading 
aa  limited  by  tha  structure* a  comllance. 

This  analysis  was  applied  by  Atkins  [47]  to  the  ease  of  axial 
splitting  of  ductile  metal  tubes.  The  ends  of  the  tubes  were  flared  by 
being  forced  down  onto  conical  mandrels,  as  shown  in  Figure  6.2.  The 
number  of  starter  cuts  produced  extremely  Interesting  axial  splitting 
behavior.  With  a  small  number  of  starter  outs,  all  cuts  would  propagate 
and  some  would  bifurcate.  Vlth  a  larger  number  of  starter  outs,  some 
cuts  would  never  develop  into  cracks.  In  all  cases,  the  final  number  of 
propagating  cracks  was  always  bs tween  8  and  12  for  50mm  diameter  and 
1.6mm  wall  thickness  mild  steel  and  aluminum  tubes.  Crack  propagation 
was  a  stable  process  with  a  reasonably  constant  radius  and  all  the 
splits  having  equal  lengths.  Atkins  applied  geometry  and  experience  to 
the  Incremental  plastic  flow  analysis  In  order  to  derive  the 
relationship! 

n  «  <8ir/  -Hr  )G  /  R)  r  t$  (6.5) 

where  a  is  the  effective  stress,  r  Is  the  radius  of  the  tube,  and  tg  Is 
the  average  hoop  strain.  Atkins’  method  predicts  n  «  12,  which 
indicates  that  the  Incremental  plastic  flow  analysis  can  be  successfully 
applied  to  some  complex  problems. 

6<5  El  MS  U  Limit  Load  Fracture  Analysis 

Crack  propagation  is  encouraged  by  transverse  strain,  and  ensuing 
necking  and  hole  nucleatlon,  in  front  of  the  crack  tip.  These  strains 
can  be  produced  by  certain  displacement  fields  within  a  fully  plastic 
structure.  The  limit  load  for  a  crack  of  given  length,  with  necking 
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neglected,  can  be  found  by  examining  the  simple  case  of  vertical  up  and 
down  cracks  extending  from  a  central  circular  hole  in  a  flat  plate. 

An  upper  bound  to  the  limit  load  is  given  by  a  displacement  field 
which  satisfies  any  displacement  boundary  conditions  and  gives  no  change 
in  volume  anywhere.  The  upper  bound  Pyb  can  be  found  from: 

$WP  -  Pub  5up  -  J  Y  «7P  «V  (6.6) 

where  Wp  is  the  plastic  work,  up  is  the  displacement  under  Pub  in  its 
direction,  Y  is  the  flow  strength,  7P  is  the  equivalent  plastic  strain, 
and  V  is  the  volume.  Experience  with  the  crushing  of  initially  deformed 
structures  suggested  separating  the  displacement  field  into  two  zones, 
as  shown  in  Figure  6.3.  The  Zone  II  displacement  field  refleots  the 
dominant  shear  strain  evident  during  tasting  which  places  the  craok 
under  tension.  Plastic  work  is  equal  to  the  applied  load  multiplied  by 
the  displacement  caused  by  plastic  deformation  of  the  structure.  the 
total  displacement  of  the  edge  experiencing  uniform  compressive  loading 
in  the  direction  of  the  applied  load  is  composed  of  ln-plane  compression 
and  drop  due  to  deflection. 

Plastic  work  is  a  function  of  stress  and  strain  throughout  the 
structure.  For  non-hardening  plasticity,  the  plastic  work  occurs  at  the 
yield  strength,  which  means  that  the  flow  stength  Y  can  be  considered 
equal  to  Since  it  is  assumed  that  the  thickness  remains  essentially 
constant  during  deformation,  Eq.  6.6  can  be  rewritten  as: 

Wp  -  h  JJ  57p  8x8y  .  (6.7) 

The  integral  must  be  evaluated  separately  for  both  zones  and  then  summed 
to  find  the  upper  bound  to  the  limit  load.  The  equations  are  simplified 
by  defining  a  displacement  ratio  of  the  maximum  vertical  deflection  to 
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the  maximum  horizontal  deflection.  Slnca  tha  aqulvalant  plastic  strains 
for  this  eaas  ars  so  complex,  the  plastic  work  is  solved  using  a 
numerical  Integration  technique, 

The  results  indicate  that  as  the  initial  flare  (dishing)  angle 
increases  the  least  uppar  bound  decreases.  Curves  of  normalized  load 
versus  displacement  ratio  for  different  initial  angles  of  flare  are 
presented  in  Figure  6.4.  The  normalized  least  upper  bound  decreases 
from  0.58  to  0,35  as  the  initial  angle  of  flare  increases  from  0  to  45 
degrees.  The  displacement  ratio  associated  with  the  normalized  least 
upper  bound  also  decreases  with  increasing  initial  angle  of  flare.' 

Hodge  [11]  and  Cordon  and  McCrum  [12]  performed  uniform 
displacement  limit  load  analyses  of  a  flat  square  plate  with  a  oentral 
circular  hole  of  varying  size,  without  any  cracks.  With  a  cutout  ratio 
(c/b)  of  40  percent,  Hodge  predicts  an  upper  bound  of  0.53  and  a  lower 
bound  of  0.41,  and  Cordon  and  McCrum  predict  an  upper  bound  of  0.60  and 
a  lower  bound  of  0.47.  The  author's  analysis  describes  a  least  upper 
bound  of  0.58  for  such  a  case.  A  structure  will  only  fail  in  the 
assumed  mode  involving  crack  propagation  if  it  does  not  first  fail  at  a 
lower  applied  load  according  to  a  different  mechanism. 

Although  the  author's  approach  is  analytically  correct,  it  does  not 
provide  any  indication  of  the  effect  of  an  initial  crack.  It  is  not 
reasonable  that  structures  will  be  unaffected  by  the  presence  of 
cracks,  regardless  of  material  properties.  It  seems  likely  that  the 
least  upper  bound  is  too  "loose,"  and  that  the  displacement  field  needs 
to  be  improved.  Further  upper  bound  limit  load  analysis  of  fractures 
could  be  conducted  using  a  displacement  field  which  specifically 
incorporates  thinning  in  front  of  the  crack  tip,  different  u  and  v 
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displacement  field*,  end  verious  hole  end  creek  alzee.  Creek 
propegetion  la  promoted  by  e  displacement  field  which  produce*  e 
combination  of  strain  and  triexlel  stress  in  front  of  the  crack.  A 


fracture  criterion  could  then  be  assumed  in  the  form  of  crack  growth  per 
unit  thinning  or  displacement  across  the  crack  tip. 

6.6  flask  EaaigifilflD  fleam  a  flitsulig  Halt 

Crack  propagation  from  a  centrally  located  circular  hole  in  a  thin 
shell  is  a  very  complex  problem.  A  J ‘based  approach  does  not  seem 
appropriate,  Atkins'  incremental  plastic  flow  analysis  can  not  be 
readily  applied,  and  a  limit  load  analysis  doss  not  provide  significant 
illumination  of  the  problem.  Further  research  is  required  in  the  field 
of  fully  plastic  crack  propagation. 
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Figure  6.1:  The  three  node*  of  loading  which  lead  to  crack  growth  [40] 


Figure  6.2:  Axial  splitting  of  ductile  metal  tube  [47] 
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Figure  6.3:  Quadrant  of  damaged  plata  and  aaaumad  diaplacamant  fialda 
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Figure  6. A.  Graph  of  upper  bound  limit  load  versus  displacement  ratio 
for  different  initial  angles  of  flare  (dishing) 


7.  EXPERIMENTS  AND  DISCUSSION 


7.1  Introduction 

An  Internal  alrblast  produces  a  variety  of  damage  inodes  examined  In 
this  paper.  The  combined  structural  effects  of  damage,  which  consist 
primarily  of  holing,  dishing  and  tearing,  have  not  been  previously 
studied.  A  small  program  of  experimentation,  consisting  of  crushing 
seven  different  specimens,  was  therefore  developed  to  gain  qualitative 
Insights  into  the  Interaction  of  various  damage  modes. 

7.2  Cumulative  Structural  Effect  of  Damage 

By  considering  the  individual  damage  effects  as  a  whole  it  is 
possible  to  paint  a  picture  of  the  residual  strength  of  structural 
plating  within  a  ship's  compartment  after  an  internal  explosion. 
Plastic  limit  load  analysis  of  a  flat  plate  indicates  that  ultimate 
strength  decreases  in  proportion  to  increasing  hole  size. 
Elastic-plastic  buckling  and  post-buckling  analysis  indicates  that 
ultimate  strength  is  little  affected  by  hole  size  until  the  holing 
reaches  half  the  plate  width,  and  is  essentially  equal  to  the  peak  load 
defined  by  the  effective  width  concept.  The  dishing  of  a  flat  plate  by 
an  explosive  loading  is  related  to  impulse  and  plate  geometry  and 
material  properties.  An  extension  of  the  concept  of  small  initial 
imperfections  to  a  dished  plats  suggests  that  the  ultimate  strength  is 
reduced  by  roughly  10  percent  from  that  of  a  flat  plate.  Finally,  the 
extensional  collapse  mode  analysis  indicates  that  the  mean  crushing  load 
of  shipboard  structural  plating  is  roughly  0.1  of  the  tensile  yield 
load.  Although  the  mean  crushing  load  decreases  with  increasing  hole 


size,  it  is  suggested  that  it  is  relatively  unaffected  by  the  extent  of 
initial  dishing.  Fully  plastic  crack  propagation  has  been  examined,  but 
its  limiting  effect  can  not  be  precisely  described.  Warship  design 
incorporates  the  use  ol  ductile  structural  plating  and  crack  arrestors 
to  negate  the  effect  of  initial  cracks  upon  residual  strength. 

Disregarding  the  potential  effect  of  crack  propagation,  the 
theoretical  combined  effect  of  damage  is  shown  in  Figure  7.1  for 
shipboard  structural  plating  with  a  length- to -thickness  ratio  of  250. 
Ultimate  strength  is  described  by  the  elastic-plastic  out-of-plane 
failure  mode  presented  in  Chapter  3.  The  plastic  in-plane  failure  mode, 
presented  in  Chapter  2,  does  not  limit  the  peak  load.  Large-scale 
dishing,  on  the  order  of  5  to  10  plate  thicknesses  from  explosive 
impulse,  only  reduces  the  ultimate  load  of  an  unperforated  plate  by  10 
percent.  Following  peak  load,  the  mean  crushing  load  ranges  from  0.06 
to  0,13  for  an  unperforated  plate,  as  presented  in  Chapter  5.  Both 
ultimate  and  mean  crushing  strength  are  little  affected  by  small  holes 
(less  than  half  the  width).  While  the  cumulative  effect  of  damage  can 
significantly  reduce  the  residual  strength  of  structural  plating,  the 
plating  retains  roughly  40  to  50  percent  of  its  original  strength.  The 
residual  strength,  even  with  extensive  damage,  is  improved  by  support 
members  which  continue  to  provide  clamped  edge  boundary  conditions  to 
structural  plating.  The  compliance  of  the  ship  tends  to  restrict 
extensive  crushing  and  provides  support  to  the  plating  to  allow  it  to 
continue  to  support  normal  service  loads. 


7.3  E&Mllamtfll  Method 

Although  not  an  exhaustive  study,  the  experimentation  was  designed 


to  examine  the  impacts  of  diffsrsnt  failure  modes  of  damaged  structure, 
end  velidete  the  enelysis  summarized  in  Chapter  7.2.  The  theoretical 
analysis  is  concerned  with  the  case  of  an  initially  flat  square  plate. 
The  simplest  way  to  perform  a  uniaxial  compressive  loading  of  the  plate 
was  to  construct  a  symmetrical  box  with  four  square  side  plates  and  two 
stiffened  end  plates.  This  arrangement  produced  a  "natural"  boundary 
condition  at  tha  side  plate  edges  which  maintained  their  right  angle 
relationship  during  Initial  ylalding,  but  did  not  prevent  out-of-plane 
deformations.  The  side  plates  were  essentially  clamped  at  top  and  bottom 
and  simply  supported  at  the  unloaded  edges. 

A  baseline  "undamaged"  specimen  was  built  of  four  undeflected  side 
plates.  Two  basic  specimen  types  ware  designed  to  reflect  different 
possible  modes  of  lateral  deflection  produced  by  tha  alrblast,  as  shown 
In  Figure  7.2.  In  order  to  model  the  different  possible  modes  of 
dishing  (lateral  deflection)  from  the  initial  undamaged  condition,  edge 
pieces  were  cut  slightly  larger  and  bent  to  shape  such  that  the  base 
dimensions  and  height  of  each  specimen  maintained  the  initial 
dimensions.  Type  A  specimens  represent  an  "Accordion  bulge"  which  has 
caused  the  edges  of  the  plate  to  deflect  along  with  the  center.  Type  B 
specimens  represent  a  "Basketball  bulga"  which  has  not  deflected  the 
edges  of  the  plate.  Due  to  geometrical  constraints  of  the  supplied 
metal,  Type  A  specimens  were  built  with  a  peak  deflection- to- thickness 
ratio  of  only  3.9,  while  the  Type  B  specimens  were  built  with  a  peak 
deflection- to- thickness  ratio  of  12.9.  Three  specimens  were  built  of 
each  type,  with  a  central  circular  hole  on  all  faces. 

The  specimens  were  built  of  1/16  inch  thick  mild  steel  (ASTM  A- 515 
grade  70)  plates  welded  together  by  the  TIG  method.  The  standard  ASTM 


character istlcs  for  this  stool  oro  presented  in  Toblo  7.1.  Tho  U.S. 
Novy  uoos  ASTM  A- 515  grodo  70  stool  in  construction  of  Novy  boilor  drums 
ond  hoodors.  During  frocturo  mechanics  tests  performed  by  tho  Navy, 
this  carbon  steel  underwent  substantial  non-linear  deformation  followed 
by  cleavage.  It  is  in  fracture  mode  transition  at  room  temperature 
[48].  The  U.S.  Navy  has  recently  adopted  HSLA-80  (ASTM  A-710  grade 
A* class  1  is  typical)  steel  for  use  as  structural  plating  in  the  hull 
and  superstructure  of  warships.  The  standard  ASTM  characteristics  for 
this  micro-alloyed,  fine-grain,  low-carbon  steel  (as  rolled  and 
precipitation  heat  treated)  are  also  presented  in  Table  7.1. 


Table  7.1:  Tensile  properties  of  ASTM  A- 515  grade  70  mild  steel 
and  ASTM-710  grade  A-class  1  low-carbon  steel  [49] 


A-515  mild  steel 

Composition  (%) 

C  <  0.27  ,  Mg  0.79-1.30 

P  <  0.035  ,  S  <  0.04 

SI  0.13-0.45 

oUT| .Ultimate  Tensile  Strength 

70-90  ksi  [485-620MPa] 

<ry  .Yield  Strength 

38  ksi  [260  MPa] 

Maximum  Thickness 

8  in  [200  mm] 

A-710  HSLA  steel 

Composition  (%) 

C  <  0.07  ,  Mg  0.40-0.70 

P  <  0.025  ,  S  <  0.025 

Si  <  0.4  ,  Ni  0. 7-1.0 

a  .Ultimate  Tensile  Strength 

90  ksi  (655  MPa) 

oy  .Yield  Strength 

80-85  ksi  (550-585  MPa) 

Maximum  Thickness 

3/4  in  (20  mm) 

Hardness  testing  (Rockwell  30T)  of  the  as -rolled  specimen  material 
yielded  Knoop  numbers  ranging  from  112  to  120  kg/mmz .  Empirical 
analysis  of  these  values  indicated  that  the  specimen  material  is  softer 


than  that  specified  by  ASTM  eharactariatlca.  Sinca  acourata  taata  of 
tha  ylald  atrangth  wara  not  conducted,  however,  tha  tanalla  yield  load 
for  tha  specimens  waa  calculated  with  tha  apaciflad  yield  atrangth 
value  of  38  ksi. 

Tha  initial  undaflactad  square  plate  had  dimanaiona  of  5.5  by  5.5 
inchaa  (140  by  140mm),  for  a  length- to* thiekneaa  ratio  of  90.  A 
warship's  atructural  plating  haa  a  length -to -thiekneaa  ratio,  neglecting 
atiffanara,  on  tha  order  of  200  to  300.  Economic  and  time  llmitatlona 
dictated  the  aalection  of  A- 515  mild  ateal  for  eonatructlng  the 
apacimena.  Although  thia  ia  a  good  ductile  ataal  to  rapraaant  the 

performance  of  HY-80  and  HSLA-80  uaad  in  U.S.  Navy  warahipa,  tha 

apacimena  were  relatively  thicker  and  aoftar  than  actual  atructural 
plating. 

A  a imp la  undamaged  apaeiman  waa  cruahad  to  aerve  aa  a  baaalina  for 
compariaon  purpoaaa.  Tha  firat  aat  of  taata  conaiatad  of  a  apaeiman  of 
each  type,  labeled  Al  and  B1  raapactively ,  with  tha  lateral  deformation 
pattern  unique  to  each.  Spacimana  Al  and  B1  had  a  2-1/4  inch  (57mm) 
central  circular  hole  on  all  aide  pletaa,  which  rapraaantad  a  cutout 
r'itlo  (c/b)  of  41  percent.  Tha  aecond  aet  of  taata  waa  aimilar  to  tha 

firat,  but  with  tha  addition  of  a  roaatte  of  8  cracka,  1/2  inch  (12.6 

mm)  long  radiating  from  each  hole ,  and  labelled  A2  and  B2 .  The  holaa  in 
apacimena  A2  and  B2  ware  tha  aame  aiza  aa  in  tha  firat  aat  of  apacimena. 
Crack  tips  wara  cut  with  hack-aaw  blades  which  had  their  aides  ground 
down  to  enaura  good  Initial  notchea  for  crack  growth.  Tha  third  aat  of 
taata  waa  designed  to  evaluate  the  effects  of  changing  hole  size  and 
initial  crack  formation.  Specimen  A3  had  tha  aame  hole  size  as  previous 
apacimena,  but  waa  modified  with  a  aat  of  short  1/4  inch  (6mm)  vertical 
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cracks  radiating  from  aach  hole.  Specimen  B3  was  modified  with  a  larger 
central  circular  hole  of  3*3/8  inches  (86mm),  which  represented  a  cutout 
of  61  percent,  but  had  no  initial  cracks.  Table  7.2  presents  a 
tabulation  of  the  initial  conditions  of  the  various  specimens.  The  six 
"damaged**  specimens  represented  a  variety  of  initial  damage  effects. 


Table  7.2:  Initial  structural  condition  for  all  specimens 


Specimen 

Baseline 

Type  A  (Accordion) 

Type  B 

(Basketball) 

un- 

damaged 

A1 

A2 

A3 

B1 

B2 

B3 

w 

Dishing  [-££) 

0 

3.9 

3.9 

3.9 

12.9 

12.9 

12.9 

Cutout  (-g-) 

0 

0.41 

0.41 

0.41 

0.41 

0.41 

0.61 

Starter  Cracks 
per  Plate 

0 

0 

8 

2 

0 

8 

fl 

The  specimens  were  crushed  on  a  Baldwin* Tate -Emery  60,000  psl 
machine  to  determine  load  versus  displacement  characteristics.  The 
testing  apparatus  is  shown  in  Figure  7.3.  Uniform  displacement  was 
regulated  by  manual  control  of  hydraulic  flow.  Load  versus  displacement 
data  beyond  ultimate  load  was  only  recorded  by  first  unloading  and  then 
reloading  slowly  until  the  specimen  resumed  deflection.  The  testing 
machine  had  a  nominal  compliance  of  1  inch  per  200,000  lb.  which  was 
considered  negligible  for  this  experimentation.  Loads  were  normalized 
with  respect  to  the  tensile  yield  load  as  described  by  Eq.  2.5: 

PY  -  A  (7.1) 


and  displacement  were  normalized  with  respect  to  the  original  specimen 
height.  Specimens  were  crushed  to  13  to  50  percent  of  their  original 


length. 


7.4  Experimental  Remit  I 

All  specimens  exhibited  generally  similar  elastic  load-deflection 
curves,  although  ultimate  strength  and  post-buckling  behavior  varied 
considerably.  The  original  load-displacement  data  for  all  specimen 
crush  tests  is  recorded  in  Apppendix  A.  All  of  the  specimens 
demonstrated  linear  elastic  behavior  during  the  initial  loading.  Loss 
of  strength  was  generally  rapid  after  the  peak  load  was  reached.  With 
continued  crushing,  however,  the  strength  of  each  specimen  leveled  off 
and  then  increased  as  local  folding  zones  collapsed. 

The  baseline  undamaged  specimen  supported  a  peak  load  of  20,600  lb. 
prior  to  collapsing,  for  a  normalized  ultimate  load  of  0.393. 
Normalized  load  versus  displacement  curves  for  all  specimens  are 
presented  in  Figures  7.4  and  7.5.  The  loads  were  normalized  with 
respect  to  the  tensile  yield  load  Py ,  and  the  deflections  were 
normalized  with  respect  to  the  original  baseline  dimensions.  The  first 
holed  specimens,  Al  and  Bl,  supported  normalized  ultimate  loads  of  0.247 
and  0.275  respectively.  The  buckling  mode  which  developed  in  specimen 
Bl  after  the  peak  load  was  reached  is  shown  in  Figure  7.6.  The  second 
set  of  holed  specimens  with  a  rosette  of  cracks,  A2  and  B2,  supported 
normalized  ultimate  loads  of  0.223  and  0.280  respectively.  Only 
specimen  A2  experienced  crack  propagation.  The  crack  growth  occurred 
only  after  the  crack  tips  had  been  re -notched  to  eliminate  the  bluntness 
which  had  developed  during  the  Initial  loading  and  after  extensive 
crushing  beyond  the  peak  load.  The  buckling  mode  which  developed  in 
Specimen  A2  after  the  peak  load  was  reached  is  shown  in  Figure  7,7.  The 


limited  creek  propagation  which  developed  wee  primarily  Mode  III  tearing 
ae  shown  in  Figure  7.8.  Specimen  A3,  with  short  vertical  cracks,  only 
supported  a  normalized  ultimata  load  of  0,253.  Although  pinching 
(necking)  was  observed  at  the  crack  tips,  specimen  A3  experienced  n a 
crack  propagation.  Specimen  B3,  with  larger  circular  holes,  supported  a 
normalized  ultimata  load  of  0.248.  The  only  visible  weld  failure 
occurred  in  specimen  23  during  crushing  beyond  40  percent  of  initial 
length.  Table  7.3  presents  a  tabulation  of  the  experimental  results. 

Table  7.3:  Ultimate  and  mean  crushing  load 
experimental  results  for  all  specimens 

Specimen  Baseline  Type  A  (Accordion)  Type  B  (Basketball) 

Un-  - 1 - 1 - 1 - 1 - 


damaged 

Al 

A2 

A3 

B1 

B2 

B3 

Ultimate  Load 
(lb.) 

20060 

12600 

11400 

12920 

14060 

14320 

12680 

Normalized 
Ultimate  Load 

0.393 

0.247 

0.223 

0.253 

0.275 

0.2B0 

0.248 

Mean  Crushing 
Load  (lb.) 

8700 

9700 

9200 

9700 

5600 

5600 

5600 

Normalized 

Mmmrt  C  v*i  i  eV  T  A 

0.17 

0.19 

0.18 

0.19 

0.11 

0.11 

0.11 

7 . 5  DigCUggjgn  af  Results 

The  experimentation  gives  some  Insight  into  the  loss  of  strength 
associated  with  damage.  The  cumulative  analysis  of  damage  effect 
failure  modes  is  applied  to  the  specimen  side  plates  and  summarized  in 
Figure  7.9.  The  experimental  results  are  also  plotted  in  Figure  7.9  for 
comparison  purposes.  Out- of -plane  buckling  caused  all  of  the  specimens 


to  collapse  at  lowar  loads  than  tha  tanaila  ylald  load.  Tha  affact  of 
damage  can,  parhapa,  ba  bast  aaan  by  comparing  damaged  specimen  strength 
characteristics  with  tha  baseline  undamaged  specimen.  Although  tha 
experimentation  consisted  of  just  seven  specimens,  some  trends  or 
observations  can  be  made  about  each  of  the  different  damage  effects. 

Boundary  conditions  affact  tha  plate  buckling  mode  and  strength 
characteristics.  Figure  7.10  shows  ho/  the  plate  buckling  coefficient 
K,  which  is  directly  ralatad  to  the  plate  critical  load,  is  affected  by 
boundary  conditions.  The  actual  boundary  conditions  imposed  upon  tha 
specimen  side  plates  is  illustrated  by  the  clamped  loaded  edges  •  C 
condition  shown  in  Figure  7.10.  Tha  plastic  limit  load  and 
elastic-plastic  analyses  presented  in  Chapters  2  and  3  only  examined 
simply  supported  edge  boundary  conditions.  The  critical  buckling  curve, 
shown  in  Figures  5.2  and  5.3,  Indicates  that  buckling  of  slender  square 
plates  of  baseline  specimen  dimensions  occurs  at  roughly  40  percent  of 
peak  load.  Although  the  actual  boundary  conditions  of  the  specimen 
plates  should  have  given  the  plates  greater  strength  than  predicted 
analytically,  out-of-plane  buckling  and  peak  load  for  tha  baseline  and 
damaged  specimens  occurred  at  lower  loads. 

All  experimental  results  for  ultimate  and  mean  crushing  load  lie 
below  the  analytic  solutions,  as  shown  in  Figure  7.9.  Since  the 
analytic  solution  for  ultimate  strength  of  an  unperforated  flat  plate  is 
known  to  be  reliable,  it  seems  likely  that  the  low  normalized  ultimate 
strength  for  the  baseline  specimen  was  produced  by  using  too  high  a 
yield  strength  value  to  determine  the  tensile  yield  load.  Hardness 
testing  indicated  that  the  specimen  material  was  softer  than  specified. 
A  lower  yield  strength  of  30  ksi  brings  the  normalized  experimental 
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results  into  11ns  with  the  analytic  solutions.  Sines  ths  true  yield 
strength  is  not  known,  however,  it  is  useful  to  compare  the  results  of 
crushing  the  damaged  specimens  with  the  baseline  specimen. 

Increased  holing  clearly  reduces  residual  strength.  The  behavior 
of  specimens  B1  and  B3  during  crushing  describes  the  effect  of  Increased 
holing.  As  the  cutout  ratio  increased  from  41  to  61  percent,  with  other 
initial  damage  effects  constant,  the  normalized  ultimate  load  decreased 
by  roughly  seven  percent  as  compared  to  the  baseline.  These 
experimentel  results  support  the  relative  loss  of  strength  with 
Increased  holing  described  by  the  elastic •plastic  analysis.  This 
indicates  that  although  holing  clearly  reduces  strength,  it  does  na£ 
cause  g  rapid  lata  fif  Strength- 

Dishing  produces  a  significant  loss  of  strength.  It  clearly  biases 
the  buckling  mode  towards  its  own  shape.  The  experimentation  Ignores  the 
effect  of  actual  thinning,  which  would  be  associated  with  dishing,  upon 
strength.  While  the  undamaged  baseline  specimen  buckled  asymmetrically 
(opposite  sides  folded  in  or  out  together) ,  all  damaged  specimens 
buckled  symmetrically  with  outward  folds  near  mid-plate.  Although  the 
Type  B  specimens  retained  greater  ultimate  strength  than  the  Type  A 
specimens,  their  post-buckling  mean  crushing  load  was  less.  It  seems 
reasonable  that  the  Type  B  specimens  retained  greater  ultimata  strength 
because  their  edges  were  Initially  straight,  while  the  Type  A  specimens 
were  initially  flared  out.  While  the  Type  A  specimens  experienced 
horizontal  folding  along  the  mid-plane  as  predicted,  the  initial 
deflection  (dished)  shape  of  Type  B  specimens  caused  the  folding  zones 
to  angle  out  from  the  central  holes  away  from  the  mid-plane.  In  the 
post-buckling  range,  dished  shape,  holing,  and  slenderness  may  become 
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the  dominant  defining  maan  crushing  load  strength. 

A  tentative  quantitative  estimate  on  the  effects  of  dishing  alone 
can  be  made  by  extrapolation  of  the  experimental  results.  Assuming  a 
linear  reduction  of  strength  with  Increased  hole  size,  the  peak  load 
results  from  crushing  of  specimens  Bl  and  B3  can  be  extrapolated  to 
determine  the  loss  of  strength  caused  by  dishing  with  no  hole.  A 
specimen  with  extensive  initial  dishing  would  thus  have  suffered  a  16 
percent  Ifigg  a£  strength  from  the  baseline  undamaged  specimen,  which  is 
a  slightly  greater  loss  than  predicted  analytically.  According  to  the 
elastic'plaatic  analysis,  however,  the  results  from  crushing  of 
specimens  Al  and  Bl  essentially  represent  the  loss  of  strength  due  to 
dishing  alone.  A  specimen  with  extensive  Initial  dishing,  but  no 
holing,  would  thus  have  suffered  a  l£i  £fi  II  percent  loss  of  strength 
from  the  baseline  undamaged  specimen.  The  effect  of  actual  thinning  of 
structural  plating  during  dishing  caused  by  an  explosion  is  considered 
negligible.  There  is  unfortunately  no  way  to  datermina  the  validity  of 
this  approach,  or  to  apply  it  to  other  dished  mode  shapes,  from  the 
limited  experimental  data. 

Both  initial  hole  size  and  deflection  shape  affect  the  mean 
crushing  load.  By  assuming  that  the  flow  strength  for  A-515  mild  steel 
equals  63  ksl,  from  Eq,  5.3  with  a  -  70  ksi,  it  is  possible  to  apply 
the  extensional  collapse  mode  analysis  to  the  test  specimens.  The  upper 
and  lower  bounds  from  the  analysis  describe  a  normalized  mean  crushing 
load  ranging  from  0.31  to  0.20,  and  from  0.24  to  0.17  for  a  specimen 
with  a  cutout  ratio  of  41  percent.  The  baseline  undamaged  specimen  had 
a  normalized  mean  crushing  load  of  roughly  0.17  during  crushing,  which 
again  suggosts  that  the  specimen  material  was  softer  than  specified. 
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The  extern tonal  collepee  mode  analysis  predicts  a  folding  zone  width  of 
1.03  Inches  for  the  baseline  specimen,  from  Eq.  5.7,  which  compares 
favorably  to  the  observed  width.  Although  the  baseline  specimen 
experienced  repeated  asymmetric  folding,  the  damaged  specimens 
experienced  only  one  major  folding  pattern  associated  with  the  initial 
dished  shape.  Type  A  and  B  specimens,  with  a  41  percent  cutout  ratio, 
experienced  different  mean  crushing  loads  of  roughly  0.19  and  0.10 
respectively.  Significant  holing  has  the  effect  of  allowing 
out-of-plane  motion  of  the  plating  around  the  hole  and  thus  reducing  the 
amount  of  energy  absorbed  by  extenslonal  deformation  of  plate  sections. 
The  extenslonal  collapse  mode  analysis  does  na£  describe  a  significant 
loss  of  mean  crushing  load  with  Increased  hole  size  until  the  cutout 
ratio  exceeds  roughly  60  percent,  which  is  supported  by  the  experimental 
results  (specimen  B3  versus  Bl).  The  initial  angled  folding  zones 
associated  with  Type  B  specimen  buckling  reduced  the  range  of  folding 
(l.e,  less  than  the  complete  folding  presumed  by  the  extenslonal 
collapse  mode  analysis) .  It  is  suggested  that  holing  and  dishing 
combined  can  reduce  the  mean  crushing  load  by  roughly  one  -  third  as  the 
dissipation  of  work  becomes  essentially  quasi -inextensional  deformation 
of  hinge  lines . 

Fracture  is  a  potentially  serious  problem  in  any  damaged 
structure,  but  was  not  a  significant  determining  factor  in  the  residual 
strength  of  the  damaged  specimens.  The  initial  starter  cuts  seemed  to 
have  little  effect  upon  strength,  possibly  because  the  fine  cuts  did  not 
have  any  necking  in  front  of  the  crack- tip  as  would  have  been  present  in 
natural  fully  plastic  fracture,  The  presence  of  starter  cracks  had  ns 
consistent  effect  upon  ultimate  strength.  The  limited  crack  propagation 
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In  specimen  A2  wee  caused  by  tearing  associated  with  complete  folding 
and  subsequent  contact  of  buckles  produced  by  extensive  crushing.  The 
mean  crushing  strength  of  the  specimens  was  little  affected  by  the 
presence  of  starter  cuts.  The  thinning  In  front  of  the  crack* tips 
indicate  that  transverse  tensile  strains  were  present,  but  insufficient 
to  promote  crack  propagation. 

A  ship's  structural  plating  is  relatively  thinner  and  harder  than 
the  specimen  side  plates,  and  may  experience  quasi-static  fracture  more 
easily  than  the  specimens  did  during  testing.  Fully  plastic  fracture  in 
complex  structures  can  not  be  described  by  current  analytic  techniques. 
Varahip  design,  however,  incorporates  the  use  ductile  materiale  to 
inhibit  fracture,  and  stringers  to  serve  as  crack  arrestors  which  limit 
any  crack  propagation. 

It  is  difficult  to  determine  the  comparative  effects  of  damage 
between  holing,  dishing  and  cracking.  The  range  of  experimentation  only 
covered  two  different  damaged  mode  shapes,  two  hole  sizes,  and  two  crack 
patterns,  as  described  in  Table  7.2.  Although  further  experimentation 
is  required  to  quantify  the  individual  and  cumulative  damage  effects,  it 
is  possible  to  apply  the  lessons  learned  from  this  research. 

7.6  Rgaidufll  S.tran&th  & £  a  Karris  After  ad  Internal  Explosion 

The  DDG-51  Arleigh  Burke  Class  destroyer  is  one  of  the  U.S.  Navy's 
most  modern  warships.  The  most  demanding  loading  any  ship 
generally  experiences  is  located  at  its  mid- section  due  to  hogging  and 
sagging  conditions  produced  by  the  ocean.  The  DDG-51  is  designed  for  a 
maximum  hogging  moment  of  176,953  ft- tons  and  a  maximum  sagging  moment 
of  87,883  ft-tons  [50].  The  maximum  compressive  loading  condition  of 
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2 

4.2  tons/in  on  the  midships  dsck  plating  occurs  during  sagging. 

From  a  structural  standpoint,  tha  worst  singla  antiship  misslla  hit 
would  ba  one  which  produces  an  internal  explosion  amidships.  An 
internal  alrblast  between  frames  220  and  254, from  a  small  antiship 
missile  such  as  the  EXOCET,  could  easily  produce  a  hole  with  a  cutout 
ratio  of  50  percent  and  dishing  on  the  order  of  10  plate  thicknesses. 
In  terms  of  longitudinal  strength,  the  deck  plating  above  the  damaged 
compartment  can  roughly  be  considered  a  10  by  10  foot  square  plate  of 
3/8  inch  thick  low-carbon  steel.  The  worst  condition  sagging  produces  a 
load  of  378,000  lb.  on  the  plating,  for  a  normalized  load  of  0.105. 
From  the  effective  width  concept,  the  undamaged  structural  plating  can 
support  a  normalized  ultimate  load  of  0.115,  neglecting  the  positive 
effect  of  support  members.  Rased  upon  the  theoretical  and  experimental 
analysis  presented  in  this  paper,  the  ultimate  strength  of  the  damaged 
structure  would  be  decreased  by  roughly  25  percent  due  to  dishing  and  5 
percent  due  to  holing,  for  a  normalized  ultimate  load  of  roughly  0.08. 
The  normalized  mean  crushing  load,  incorporating  the  effect  of  holing 
but  neglecting  dishing,  would  range  from  roughly  0.07  to  0.04.  These 
strength  characteristics  indicate  that  the  plating  can  not  support 
extreme  sagging  conditions  prior  to  ultimate  load  on  its  own,  but  could 
probably  continue  to  support  normal  service  loads  before  and  after  peak 
load.  A  local  failure  of  the  structural  plating  in  one  or  more 
compartments  does  not  mean  that  the  ship  will  experience  global  failure. 
The  ship's  support  members  and  structural  compliance  would  offer 
continued  support  to  the  plating  to  prevent  local  collapse. 


97 


.4  - 
-O  - 
0J2  - 
OJ  - 


von  kArmAn  effective  width  for 
UNPERFORATED  PLATE 


LOSS  OF  STRENGTH 
DUE  TO  DISHING 


ELASTIC »PUSTIC  ULTIMATE  LOAD 
OUT -OF- PLANE 


MEAN  CRUSHING  ^2Z222Z2Z2Z2ZZ1 


— T— 
0£ 


0.4 


06 

CUTOUT  RATIO  [-£-) 


Figure  7.1:  Summary  of  damage  effect  failure  modes  for  shipboard 

structural  plating  with  a  length- to- thickness  ratio  of  250 
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Figure  7.2:  Specimen  types  with  different  simplified  dishing  (lateral 
deflection)  modes 
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Figure  7.3:  Photograph  of  tenting  apparatus 
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Figure  7.4:  Graph  of  normalized  load  versus  deflection  for  baseline  and 
Type  A  (Accordion  bulge)  specimens 
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Figure  7.7:  Photograph  of  specimen  A2  shortly  aftar  ultimata  load  was 
raachad 


Figure  7.8:  Photograph  of  tearing  from  a  starter  crack  in  specimen  A2 


Figure  7.9:  Summary  of  damage  affact  fallura  mod**  and  axparlmantal 
raaults  for  Gapaciman  alda  plataa  with  a  length- to- thickneaa  ratio  of  90 


Figure  7.10:  Graph  of  plata  buckling  coefficient  versua  aapact  ratio 
for  different  boundary  condition*  [18] 
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8 .  CONCLUSION 


Shipboard  atructural  plating  can  continua  to  aupport  normal  service 
loading  even  after  extensive  damage.  An  internal  explosion  produces 
tearing,  holing,  and  dishing  of  structural  plating.  The  shipboard 
structure  and  damage  effects  were  simplified  in  this  paper  in  order  to 
analyze  the  individual  damage  effects  and  determine  how  they  Interact. 
A  small  program  of  experimentation  roughly  validated  the  cumulative 
analysis  of  damage  effect  failure  modes. 

Structural  collapsa  is  a  function  of  mechanical  properties, 
structural  geometry,  and  collapse  mode.  (1)  The  ultimate  strength  of 
structural  plating  is  strongly  affected  by  slenderness.  (2)  The 
elastic-plastic  analysis  accurately  describes  the  ultimate  load  for 
slender  plates,  and  indicates  that  small  holes  (cutout  ratio  less  than 
50  percent)  have  little  affect.  (3)  Dishing  on  the  order  of  5  to  15 
plate  thicknesses  reduces  ultimate  strength  by  roughly  15  to  35  percent. 
(4)  The  buckling  mode  involves  folding  as  decrlbed  by  the  extenslonal 
collapse  mode  analysis.  (5)  Tha  normalized  mean  crushing  load,  which 
describes  the  average  compressive  strength  maintained  after  the  ultimate 
load  has  been  reached,  is  roughly  0.1  for  undamaged  shipboard  structural 
plating,  and  ranges  from  0.05  to  0.1  for  significantly  damaged  plating. 
(6)  Although  crack  propagation  was  limited  during  crushing  of  the 
damaged  specimens,  fully  plastic  fracture  could  be  more  significant  in 
shipboard  structural  plating.  Continued  reserach  is  required  to  further 
quantify  these  conclusions  and  apply  the  results  to  Improved  ship 
designs . 


8.1  Ideal  for  Continued  Research 

This  paper  has  presented  the  first  qualitative  theoretical  and 
empirical  analysis  of  the  cumulative  effect  of  damage  upon  structural 
plating.  There  are  many  avenues  of  continued  research  required  to 
develop  a  comprehensive  quantitative  understanding  of  damage  effects. 
Ideas  for  further  research  Include:  (1)  determine  the  effect  of  ship's 
compliance  upon  plate  boundary  conditions,  (2)  determine  the  effect  of 
stiffeners  upon  various  damage  modes,  (3)  Investigate  rectangular  plate 
geometries,  (4)  determine  the  deflection  of  supporting  members 
associated  with  plate  geometries,  (5)  examine  dynamic  effects,  (6) 
examine  fully  plastic  fracture,  (7)  conduct  extensive  critical 
experimentation  to  determine  quantitative  relationship  between  damage 
effects.  It  is  important  that  any  testing  properly  model  shipboard 
structural  geometry  (such  as  b/h)  and  material  characteristics  (such  as 
0y/E) .  The  author  suggests  that  determining  the  effete  of  fully  plastic 
fracture  and  shipboard  structural  compliance  are  particularly  Important. 
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